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INTRODUCTION 



The J.H. Campbell Plant is comprised of three coal- 
fired operational units. Units 1 and 2 use Lake Michigan 
water drawn through Pigeon Lake for cooling purposes. Unit 
3 draws cooling water from an intake located in 
approximately 1 1 m of water (1.1 km offshore) in Lake 
Michigan. Beginning in September 1980, heated water from 
all three units was discharged through a common offshore 
structure located in approximately 6 m of water (0.3 km 
offshore) . 

Concurrent with reports on larval, juvenile, and adult 
fish (Jude et al. 1978, 1979, 1980, 1981, 1982), this is the 
last in a series of five reports on the benthos that began 
in 1977. Previous reports (Jude et al. 1978; Winnell and 
Jude 1979, 1980, 1981) have dealt with annual, monthly, 
depth, and regional density fluctuations of the benthos and 
grain size analysis of substrates occurring in the vicinity 
of the Campbell Plant from the perspective of each 
particular year. Comparisons were made to determine the 
relevance and longevity of regional trends as related to 
plant operation and general lake-wide distribution (for more 
on lake-wide distributions in addition to Winnell and Jude, 
see Powers and Robertson 1965; Robertson and Alley 1966; 
Hiltunen 1967; Alley 1968; Mozley and Garcia 1972; Mozley 
and Alley 1973; Mozley 1974, 1975; Alley and Mozley 1975; 
Mozley and Winnell 1975; Mozley and Howmiller 1977; Nalepa 
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and Quigley 1980; Nalepa and Robertson 1981). When compared 
with lake-wide distributions, we have observed some 
similarities (dominance of amphipods, oligochaetes, and 
chironomids) and dissimilarities (very high abundance of 
Pontoporeia hoyi ) . Other investigations near the Campbell 
Plant have indicated some regional differences (Truchan 
1970; Consumers Power Company 1975). We noted during 
preoperational years (1978-1980) that, while depth and 
temporal factors accounted for much of the variability in 
benthic distributions, regional differences did exist and 
needed to be monitored. With the inclusion of the first and 
only full operational year for which benthic data will be 
collected, we were able to employ a modification of the 
statistical design proposed by Winnell and Jude (1980). The 
analysis of variance (ANOVA) design was modelled after a 
similar 8-yr design (Johnston 1974) used to test for a plant 
effect at the D.C. Cook Nuclear Power Plant located 113 km 
south along the shoreline of Lake Michigan [unpublished 
data, Great Lakes Research Division (GLRD)]. The modified 
statistical design was applied to benthic density 
fluctuations found in the 3 to 15-m depth range sampled from 
1978-1981 near the J.H. Campbell Power Plant. We found no 
differences among density fluctuations at the ANOVAs' levels 
of sensitivity. In addition, we concluded that entrainment 
of malacostracans from pumping of cooling water via intake 
structures located in Lake Michigan had no apparent effect 



on lake populations of malacostracans or on overall lake 
ecology. 

METHODS 

LAKE MICHIGAN BENTHIC FIELD SURVEYS 
Benthos Field Survey Design 

The survey was composed of 10 stations located along 
two transects perpendicular to the eastern shoreline of Lake 
Michigan near the J.H. Campbell Power Plant, Ottawa County, 
Michigan (Fig, 1). Along each transect, stations were 
located at 3-, 6-, 9-, 12-, and 15-m depths. The first 
transect represented the treatmient area (inner region) 
located 0.16 km north of the offshore intake and discharge 
structures. The second transect represented the reference 
area (outer region) located 5.0 km north of the offshore 
intake and discharge structures. The survey design was the 
same during 1979, 1980, and 1981. 

Benthos Collection and Processing 

Benthic macroinvertebrate samples were collected on 13 
April, 13 July, and 8 October 1981. Sixty samples per month 
were collected from the University of Michigan's R/V Mysis. 
Sediment samples, which were collected during 1979 and 1980, 
were not collected during the 1981 field season. 

Benthic macroinvertebrates were collected with a 
triplex (three-chambered) Ponar grab sampler (Mozley and 
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Chapelsky 1973). Each chamber of the Ponar samples 0.0165 

2 

m . Contents of one side chamber were used to estimate 

numbers of benthic macroinvertebrates occurring in a square 
meter of lake bottom. Six replicates were collected to 
estimate benthic populations at a station on each sampling 
date. 

The portion of the sample used to estimate benthic 
macroinvertebrates was placed in a "funnel-shaped hopper" 
(see Mozley 1975 for details) aboard the ship. Samples were 
washed through a 0.2-mm mesh net to concentrate animals and 
remove excess sediment and debris. Concentrated samples 
were stored in labelled 0.5-liter Mason jars and preserved 
with carbonate-buffered, 4% formaldehyde solution. Samples 
were returned to the GLRD Benthos Laboratory for sorting and 
identification. 

Sorting and initial identification of organisms were 
performed using dissecting microscopes (3-30X). Specimens 
unidentified at the genus/species level were mounted on 
slides with Amman's lactophenol clearing medium and 
identified using compound microscopes (40- 1 , OOOX) . 

Initial generic identification of chironomids was 
determined using an unpublished trial key to the chironomids 
(A.L. Hamilton and O.A. Saether, personal communication, 
Freshwater Institute, Winnipeg, Manitoba, Canada and 
Zoological Museum and Department of Morphology, Systematics 
and Animal Ecology, University of Bergen, Bergen, Norway). 
Where species were determined for chironomid genera. 



"cf." refers to uncertain larval identification at the 
species level. Most species designations concur with reared 
specimens from the D.C. Cook Nuclear Power Plant, 
southeastern Lake Michigan (see Mozley 1975), which are 
maintained in the GLRD Benthos Laboratory's permanent 
collection. However, as none of the chironomid larvae from 
the Campbell Plant have been reared, identifications at the 
species level have been assigned the uncertainty designator 
"cf.". The designator "gr." refers to a "group" of species 
undeterminable from larvae and was associated with the 
genera Chironomus and Paracladopelma . Morphology and 
taxonomy of other chironomid genera and species were 
determined from the following references: Lenz (1954), 
Roback (1957), Curry (1958), Beck and Beck (1969), Saether 
(1969, 1971, 1973, 1975, 1976, and 1977), Hirvenoja (1973), 
Maschwitz (1975), Jackson (1977), and Soponis (1977). 

Naidids were identified using Hiltunen's key to the 
naidids (see Hiltunen and Klemm 1980). Tubificids were 
identified using an unpublished key to aquatic oligochaetes 
of the Great Lakes (J.K. Hiltunen, personal communication. 
Great Lakes Fishery Laboratory, U.S. Fish and Wildlife 
Service, Ann Arbor, Michigan). Gastropods and pelecypods 
were identified using a key to the mollusks of the Great 
Lakes (Mackie et al. 1980). Amphipods were identified using 
the keys developed by Holsinger (1976) and Pennak (1978). 



statistical Design, Sensitivity of the ANOVA , 
and Determination of a Plant Effect 

The statistical design initially proposed to determine 
plant effect (in terms of heated effluent discharged into 
Lake Michigan) and relative sensitivity of the design to 
population changes was a balanced mixed-model, nested 
analysis of variance (ANOVA). This design originated from 
an 8-yr design developed by Johnston (1973, 1974) and 
employed at the D.C. Cook Nuclear Power Plant (unpublished 
data, GLRD). The present application was based on a 4-yr 
design with 2 preoperational yr (1978 and 1979) and 2 
operational yr (1980 and 1981) anticipated, making our 
design a balanced modification of the 8-yr design. A 
balanced design was necessary not only for computational 
purposes, but for more reliable conclusions. However, as 
Unit 3 at the Campbell Plant did not become operational 
until September 1980, it was necessary to alter the 
statistical design due to the resulting unbalanced nature [3 
preoperational yr (1978-1980) and only 1 operational yr 
(1981)]. This complication was previously noted by Winnell 
and Jude (1981). The analysis of variance for each taxon 
was performed on log^Q(x+1) transformed densities (Elliot 
1971) by the University of Michigan's AMDAHL 470V/8 computer 
using BMD8V, a program of the BMD Statistical Software 
facility supported by the Statistical Research Laboratory at 
the University of Michigan. 

The original design was structured such that there were 
five factors: Construction Time [before and after 
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discharging of heated effluent via offshore discharge 

structures began], Year [1978-1981], Region [treatment 

(inner) and reference (outer)], month [April, July, and 

October], and Depth [3 m, 6 m, 9 m, 12 m, and 15 m] . All 

factors were fixed factors, except Year, which was random 

and nested within Construction Time. Because the 

Construction Time factor was no longer a balanced factor, 

thereby complicating computation and interpretation of 

results, it was necessary to eliminate this factor and 

utilize the Year factor as an unnested factor. In order to 

do this, the design of the ANOVA was altered to a four-way 

mixed, factorial model where all factors were fixed, except 

Year, which remained random. 

While determination of a plant effect was measured by 

the F-ratio of the mean square errors of the Construction 

Time x Region interaction and Year x Region interaction 

(MSpp/MSyp) in the original design, elimination of the 

Construction Time factor required utilization of a different 

test for plant effect. The test for plant effect that we 

employed in the new design was directly related to the 

degree of sensitivity of the ANOVA. The ANOVA' s sensitivity 

to population changes was determined by R, the least 

detectable true ratio (see Cohen 1969; Johnston 1974; 

Winnell and Jude 1980): 

^0.5. -oO'Sr 

R > 10"^ ^ and R < 10 ^ ^ 



The quantity 6 (least detectable true difference) in the R 
equation was determined by Johnston (1974) by re-expressing 
the formula of Sokal and Rohlf (1969): 

6= s(2/n)°-5 (t^[^j * t2,,.p)[^3)2, 

where; 

6 = least detectable true difference, 
s = true standard deviation, 
Y = degrees of freedom, 
n = number of observations, 
t = Student' s t , 
a = significance level , 

P = desired probability that a difference 
will be found to be significant. 

Calculations of 6 from Campbell ANOVAs assumed: 



s = mean square error, 

Y = infinity, 

n = average number of observations per 
cell in the reduced ANOVA [the total 
number of observations/4 (see below)], 

a = 0.05, 

P = 0.95. 



Reducing the number of ANOVA cells such that columns were 
years pooled into preoperational (1978-1980) and operational 
(1981) years, and rows were regions (inner, outer), there 
remained four cell means corresponding to those of Johnston 
(1974, p. 27). 

Subsequent substitution of 6 into the R equation 
permitted us to determine the minimum amount of change 
necessary to detect that a change had in fact occurred. 
This change corresponded to the sensitivity of the ANOVA, 
calculated at the 5% level of significance and with 95% 
certainty that a difference at least equivalent to 6 would 
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be regarded as significant. Consequently, the least 
detectable true ratio expressed the minimum degree of change 
(either increase or decrease) which the inner region 
population density for a given animal must undergo relative 
to a similar estimate from the outer region in order to 
detect plant effect using the four-way mixed, factorial 
ANOVA. Therefore, the quantity R not only described the 
sensitivity of the ANOVA, but delimited the threshold below 
which, we can conclude, there was either no plant effect or 
the effect of plant operation did not exceed the ANOVA' s 
level of sensitivity (see RESULTS AND DISCUSSION for 
detail) . 

To determine if the respective R values were exceeded 
by density changes observed at Campbell from 1978-1981, the 
four cell means from the reduced ANOVA were contrasted 
according to Johnston (1974): 

R' = (IA+1/0A+1)/(IB+1/0B+1), 
where; 

R' = actual abundance change ratio, 

-2 

lA = true mean number of animals m 

(geometric mean) in the inner region 
after plant began operation, 

-2 

OA = true mean number of animals m 

(geometric mean) in the outer region 
after plant began operation, 

-2 

IB = true mean number of animals m 

(geometric mean) in the inner region 
before plant began operation, 

-2 

OB = true mean number of animals m 

(geometric mean) in the outer region 
before plant began operation, 
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Therefore, our test for plant effect on benthic populations 
as reflected in density fluctuations compared R' with R. If 
R' £ 10 "^ ^ and R' > 10^ ^ , we were 95% certain the 
observed degree of change was significant at the 5% level, 
i,e., for a given animal, the inner region density relative 
to the outer region density significantly increased or 
decreased. This condition was judged to strongly indicate 
an effect due to discharge of heated effluent from plant 
operation. 

Conversely, if R' > 10 ^ ^ and R' < 10^ ^, we were 
95% certain the observed degree of change was not 
significant at the 5% level of significance, i.e., for a 
given animal, the inner region density relative to the outer 
region density did not experience a significant increase or 
decrease in abundance. This condition could be judged as 
either no effect due to discharged heated effluent, or a 
non-detectable effect (see RESULTS AND DISCUSSION for 
detail) . 



ARTIFICIAL SUBSTRATE STUDY ON THE RIPRAP AREA 

The riprap, comprised of large, irregular pieces of 
crushed limestone (0.5 to 2.0-m diameter), covers the intake 
and discharge structures. This protective cover extends 
from the shoreline to a depth of 11 m, ranges in width from 
9-18 m, and has a total estimated area of 52,400 m^ (Jude et 
al. 1982). Artificial substrates used to estimate both 
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benthic diversity and density on the riprap followed those 
utilized by Lauritsen (1979). As the geometric design of 
the artificial substrates was a truncated cone, area of 
individual substrates was calculated by: 

A = p[r^ + R^ + (rR)s3 , 

where; 

A = area of a truncated cone, 

p = 3.1416, 

r = radius of circle at 

top of cone, 
R = radius of circle at 

bottom of cone, 
s = slant height of cone. 

Each artificial substrate was constructed of quick drying 
concrete formed in paper cups, such that when dried, r = 

26 mm, R = 36 mm, and s = 88 mm, with total area of each 

2 

individual artificial substrate determined to be 0.02334 m . 

Five artificial substrates were placed in each of five 
wire baskets. The contents of each basket were considered a 

replicate. The total area sampled per wire basket was 

2 2 . 

equivalent to 5 times 0.02334 m or 0.11668 m , assuming 

that 100% of each artificial substrate surface area was 
available to potential benthic colonizers. However, as 
positioning of the five artificial substrates relative to 
each other did not ensure that all available surface area 
was exposed to colonization, we followed the methodology of 
Lauritsen (1979) by assuming 10% of the area was unavailable 
to colonization due to posturing of substrate surface areas 
against themselves (0.11668 m^ - 0.01167 m^ = 0.010501 
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2 

m ). Therefore, the conversion factor used to convert the 

-2 

number of individuals per wire basket to number m was 

9.52. 

Wire baskets containing the artificial substrates were 
lodged by divers in the riprap surrounding the intake 
structure on 15 July and were retrieved on 22 October. 
Scuba divers placed dislodged wire baskets immediately in 
0.1-mm mesh bags, which were sealed and returned to the 
ship. On board ship, each bag and its contents were placed 
in separate 18.9~liter plastic buckets, preserved with 4% 
formaldehyde solution, and capped. Samples were returned to 
the GLRD Benthos Laboratory for sorting and identification 
procedures following those previously described for 
processing benthic samples collected by Ponar grabs. 

MALACOSTRACANS ENTRAINED AND COLLECTED BY NET AND SLED TOWS 

Malacostracans collected in net tows were sampled using 
a 0.5-m diameter, nylon plankton net of no. 2 mesh (363- 
micron aperture). A Rigosha flowmeter attached to the 
center opening of the plankton net was used to calculate 
volume of water sampled. When flowmeters were not available 
or stopped functioning, average flowmeter values were 
computed from readings available from the same stations at 
other times or from stations of comparable depth. When 
flowmeter readings were conspicuously different from other 
tows at the same station, averages of readings for the 
appropriate station and diel period were used. All meter 
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revolutions were converted to volume filtered using 1 
revolution = 15 liters of water. Flowmeters were calibrated 
in a swimming pool by walking a measured distance with a 
flowmeter attached to a 0,5-m diameter hoop without the net 
(see Jude et al. 1979) . 

Net tow samples from the beach were collected at 
station R (in the vicinity of the intakes) in Lake Michigan 
(Table 1). Three people simultaneously hand-towed two nets 
for a distance of approximately 61 m once during the day and 
once at night. These tows were performed once per month 
during April and September and twice per month in May, June, 
July, and August during 1981. 

At deeper stations, malacostracans were collected in 
horizontal, 5-min tows taken at discrete depths parallel to 
shore at six stations (I, J, L, N, 0, W) in Lake Michigan 
(Table 1). Sampling was conducted during the day and night 
on the same schedule as beach samples. A total of 476 net 
tow samples were collected during 1981, with only the 7 May 
station R sample missing (0.5-m night net tow). 

Net tows at depths less than 3 m were taken from 

outboard motorboats. The University of Michigan's R/V Mysis 

was used for tows at deeper stations. For each tow, the 

following procedure was implemented: 

1) Plankton net (with attached 0.47-liter 
Mason jar) and depressor were lowered 
to desired depth (average ship speed 
was 3 to 6 km/h) . 
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Table 1. Net and sled tow stations and respective 
sampling depths (m) at which malacostracans were 
collected during 1981 near the J.H. Campbell Plant, 
eastern Lake Michigan. The last row in the table 
denotes bottom depth at each station. Sled tow samples 
were collected on the bottom, while net tow samples 
were collected at all depths noted in the table, except 
at the bottom depth. 

Station 
R I J L N W 



0.5 0.5 



1.0 1.5 



0.5 


0.5 


0.5 


0.5 


0.5 


2.5 


2.0 


2.5 


3.0 


4.5 




4.0 


4.5 


6.0 


8.5 




5.5 


6.5 


9.0 


11.5 






8.5 


11.0 


14.0 


3.0 


6.0 


9.0 


12.0 


15.0 
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2) Plankton net was towed horizontally 
for 5 min when the desired depth was 
was attained. The amount of cable or 
rope needed to obtain the desired 
depth was calculated trigonometr ically . 

3) Plankton net was hauled to surface where 
contents were washed into the attached Mason 
jar using a water hose. 

4) Contents of the Mason jar were preserved 
with 40 mL of buffered formaldehyde, 
labelled, and sealed. 



Sled tows along the lake bottom were made with a 
benthic fish larvae sled equipped with a flowmeter (Yocum 
and Tesar 1980). A single 5~min sled tow was performed once 
during the day and once during the night at all Lake 
Michigan stations coincident with net tows as time and 
weather permitted. All sled tows were made from outboard 
motorboats. A total of 140 sled tow samples were collected 
during 1981. • 

Due to the design of the Unit 3 cooling system, 
entrainment samples were collected at the opening of the 
5.5-m (18-ft) diameter intake pipe which enters the 
uncovered intake canal near the Lake Michigan shoreline. 
Collection at this point ensured that very few, if any, 
malacostracans occurring in the intake canal were sampled, 
which could potentially bias entrainment samples with 
malacostracans not originating from Lake Michigan. 

Entrainment samples were collected from a raft situated 
directly above the juncture where the buried intake pipe 
carrying cooling water from Lake Michigan entered the on- 
land intake canal. A O.S-m diameter plankton net, identical 
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to that used for field sampling {363-micron mesh), was 
lowered into a central position in the mouth of the intake 
pipe. An IS-kg weight was needed to keep the net at the 
desired depth. The amount of cable necessary to maintain 
this central position was calculated trigonometrically . 
Location of the raft and amount of cable necessary to ensure 
correct positioning of the net in the opening was verified 
by scuba diver observations. 

Four replicate samples were collected four times during 
each 24-h period sampled (day, dawn, dusk, and night). Each 
replicate consisted of a 10~min sample. Sampling was 
performed on a weekly basis from May through September, 
three times per month in April and October, and twice per 
month from January through March and November through 
December 1981. The periods of dawn and dusk were defined as 
the period extending from 1 h before to 1 h after both 
sunrise and sunset as determined from the Nautical Almanac 
Office (1981), United States Naval Observatory, Washington 
D.C. During 1981, a total of 590 entrainment samples was 
collected, with only the four entrainment samples from 10 
July missing. However, as the accuracy of entrainment 
sample flowmeter readings of <500 revolutions (<7.5 m 
filtered) was suspected to be low due to fouling or 
malfunction, 48 of the 590 entrainment samples collected 
were considered unreliable and were not included in the 
dataset. 
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All samples were returned to the GLRD Fishery 
Laboratory where malacostracans (Amphipoda, Mysidacea, and 
Isopoda) were removed from each sample and stored in 4- or 
8-mL vials. For samples having excessively large numbers of 
amphipods or mysids, the first 50 individuals of each taxon 
were removed, with the remainder enumerated to provide the 
total number in the sample from which an appropriate 
conversion factor was derived. At this stage in processing 
malacostracans, samples were entrusted to the Benthos 
Laboratory for identification and measurement. 

Malacostracans were identified to the lowest practical 
taxonomic level. For each specimen, body length was 
measured to the nearest 0.1 mm: for amphipods, from the 
insertion point of the first antenna on the head to the tip 
of the telson; for mysids, from the tip of the rostrum to 
the posterior extension of the telson; and for asellids, 
from the anterior portion of the head to the tip of the 
uropods. In addition, head length was determined for 
numerous, but not all, specimens. Reproductive status 
determined for all specimens was straightforward (juvenile, 
gravid female, spent female, and male), except for Mysis 
relicta , which followed the method of McWilliam (1970). 

Malacostracans were sorted, identified, assigned a 
reproductive status, measured, and grouped into 0.5-mm size 
categories. These data along with appropriate station 
information (e.g., depth, time, water temperature), were 
entered directly onto an internal computer datafile of the 
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Michigan Interactive Data Analysis System (MIDAS) at the 
University of Michigan Computing Center. With the aid of 
data analysis routines of MIDAS, appropriate statistics were 
generated for malacostracans occurring in net tows, sled 
tows, and entrainment samples. 

Because net tow sampling involved sequential vertical 
subsamples of depth at each station, density and meter 
revolutions for each subsample were corrected to account for 
any contamination from individuals occurring at other depths 
through which the net had to be hauled when retrieved and to 
provide a more accurate measure of the number of cubic 
meters of water passing through the net. The adjustment 
procedure for organisms captured in all tows other than 
surface tows is outlined in detail in Figure 6 of Jude et 
al. (1982). The method consisted of sequential subtraction 
of numbers of malacostracans from the lower water depth 
levels based upon densities observed in upper water strata. 
It was assumed that malacostracans were homogeneously 
distributed within a water stratum and that nets passing 
through a particular stratum from a lower level would catch 
individuals in proportion to the volume of water filtered. 
Organisms from all tows conducted below the surface stratum, 
which were probably caught during the vertical haul 
following termination of the horizontal tow, were removed 
via calculation from the final total malacostracan 
densities. It was assumed that any contamination incurred 
while lowering the net was negligible. The effect of 
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differential vertical distribution due to size was mitigated 
by stratifying malacostracans from each sample into 0.5-mm 
length intervals. 

Biomass estimates for Pontoporeia hoyi and Gammarus 
spp. were generated by summing the products of the weights 
for mid-point lengths for each successive 0.5-mm size class 
times the number of individuals occurring in each respective 
size class. The length-weight curve (ash-free dry weight in 
mg) used for P. hoyi and Gammarus spp. was developed by 
Johnson and Brinkhurst (1971) for Lake Ontario specimens: 

P. hoyi wt(mg) = 0.014L^-^^°, 

2 444 
Gammarus fasciatus wt(mg) = 0.014L 

Similar biomass estimates determined for Mysis relicta were 
made utilizing the length-weight curve (dry weight) of 
Morgan (1976): 

Mysis relicta wt(mg) = 0. 00016L^ * "^^ . 

Biomass was not estimated for remaining taxa. 

As previously noted, pumping operations at Unit 3 draw 
cooling water from a depth of 1 1 m ( 1 . 1 km offshore) in Lake 
Michigan. In the process, many types of animals are 
entrained, many of which are malacostracans. In an attempt 
to evaluate the effect of entrainment on populations of 
malacostracans in the lake, we assumed the area of greatest 
potential impact was a 2-km area encompassing the 3 to 18-m 
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depth regime. This area is equivalent to 4 x 10^ m^. The 

volume of water (V) overlying the area with these dimensions 

is 4.2 X lo'^ m"^ [V. ^ T = V ^ , + V . = (3 m)(2,000 

total rectangle prism xu/v^^uuu 

2 2 

m) + 0.5(15 m) (2,000 m) ]. Both area and volume estimates 

were assumed to be constant during any given time period. 
In addition to these, we calculated an average daily pumping 
rate for Unit 3 based on data supplied by Campbell Plant 
personnel for several different time periods of interest. 
The maximum daily pumping rate for Unit 3 is 1.766 x 10^ m'^ 
day f which we used to estimate maximal entrainment. 
However, the pumping rate which we assumed was the average 

actual daily pumping rate during any given day during 1981 

6 3 — 1 
was 1.115 X 10 m day (based on pumping data from January 

through December) and was used to estimate actual 

entrainment. Finally, when utilizing net tow densities for 

which estimates were limited to the April through September 

time period, we used net tow density averages from the 3 to 

15-m stations only and the actual average daily pumping rate 

6 3 — 1 
of 1.177 X 10 m day . Using these assumptions in 

conjunction with measured densities at specified time 

periods, we were able to evaluate the effect of entrainment 

on lake abundances of certain malacostracans using the 

following equations: 

1) Total average number of P. hoyj or 
M. relicta present in the 2~km area 
in the 3 to 18-m depth regime at any average, 
given time during the year 
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p. hoyi = [4 X 10^ m^][3,849 m ^ 

(Inner region density at 

3 to 15 m, 1981)] 

1 3 
= 1 .54 X 10 '^ m-^. 



M. relicta = [4 x 10^ m^][188 m"^ 



(Morgan and Beeton 1978)] 
= 7.52 X 10^. 



2) Annual entrainment based on January-December 
average entrainment for each taxon 

Maximum number entrained annually = (no. entr. m ) 

(1.766 X 10^ 
m days ) 
(365 day yr'^ 



3 
Actual number entrained annually = (no. entr. m ) 

(1 .115 X 10^ m^ 
m"^ day"'') 
(365 day yr""''). 



3) The percentage that the density of each taxon 
entrained represents of its respective benthic 
density in the 2-km area 

Maximum % = (max. no. entrained annually/ 
total no. in 2-km area) (100). 



Actual % = (actual no. entrained annually/ 
total no. in 2-km area) (100). 
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4) The percentage that the density of P. hoyi 

and M. relicta (April-September) represents of their 
respective average daily pelagic densities in the 
volume of water overlying the 2-km area 



Maximum daily % (P. hoyi ) = 

[(1.766 X 10^ m^ day"'')(2.347 m"^)/ 
(4.2 X lo"^ m^)(1.098 m"-^)](100) 
= 8.99%. 

Maximum daily % (M. relicta ) = 

[(1.766 x 10^ m^ day~'')(0.1012 m~^)/ 
(4.2 x lo"^ m^ day""*) (0.0569 m"^)](100) 
= 7.48%. 



Actual daily % (P. hoyi ) = 

[{1.177 X 10^ m^ day"'')(2.347 m~^)/ 
(4.2 X 10^ m-^ day"''){1.098 m"^)](100) 
=5.99%. 



Actual daily % (M. relicta = 

[(1.177 X 10^ m^ day"''){0.1012 m"^)/ 
(4.2 X lo'' m^ day""*) (0.0569 m"^)](100) 
= 4.98%. 
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5) The percentage that densities of pelagic P. hoyi 
represent of densities of benthic P. hoyi in the 
2-km area at any average time during April 
through September 



% pelagic = [(4.2 x 10^ m'^)(1.098 m~^)/ 
1 .54 X 10^°](100) 
= 0.299%. 



Maximum and actual biomass of malacostracans entrained 
on an annual basis was calculated in the same manner as was 
done for malacostracans entrained annually, except the 

average number of mg m was substituted into the equation 

2 

in place of no. m . For P. hoyi , this quantity was 1.272 mg 

-3 . -3 

m ; for mysids, 0.0296 mg m ; and for Gammarus spp. , 

-3 
0.0788 mg m . In the case of P. hoyi , accurate benthic 

abundances for individuals occurring in the size categories 

< 3 mm, 3 to 5 mm, 5 to 7 mm, and >7 mm were known from 

Ponar grabs collected at 3 to 15 m in the inner region 

during April, July, and October 1981. Based on the 1981 

inner region benthic abundance of P. hoyi in each size 

category and the weight (calculated from Johnston and 

Brinkhurst 1971) of an individual at the mid-point of each 

size class (i.e., 2 mm, 4 mm, 6 mm, and 8 mm), average 

-2 

benthic biomass was 1,541.9 mg m . Therefore, the biomass 

of the benthic P. hoyi population at any average, given time 

6 2 
during the year was equivalent to (4 x 10 m ) (1,541.9 mg 

-P Q 

m ) = 6.168 X 10 mg = 6168 kg. Therefore, the percentage 
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that the biomass of entrained P, hoyi comprised of the 
benthic biomass is given by: 

Maximum annual % = (Maximum annual biomass 
of entrained P. hoyi /6, 168 kg)(100). 

Actual annual % = (Actual annual biomass 
of entrained P. hoyi /6, 168 kg) (100). 

Estimating the effect of entrainment on Gammarus spp. 
was limited to October as benthic abundance on the riprap 
was representative of October only. The product of the 

average Gammarus spp. abundance from artificial substrates 

-2 2 

(139 m ) and the area of the riprap (52,400 m ) represented 

the estimated total benthic Gammarus spp. density available 
to entrainment (7.284 x 10 individuals). As the October 

entrainment density of Gammarus spp. at the 3- to 15-m 

~3 
stations was 0.2005 m , percent entrained was determined 

such that; 



Maximum daily % entrained = 
[(1.766 X 10^ m^ day"'')(0.2005 m~^)/7.284 x 10^](100) 

= 4.86%. 

Actual daily % entrained = 
[(0.948 x 10^ m^ day" b (0 .2005 m"^)/7.284 x 10^](100) 

= 2.61%. 



25 



The above constitute the general, and in several cases, 
the specific equations used to estimate the impact of 
entrainment on malacostracans. Quantities generated from 
these equations appear in the section entitled "Impact of 
entrainment on malacostracans." 

Finally, it should be noted that annual entrainment 
estimates from the Cook Plant have been calculated with the 
same methodology used at Campbell, which was based on 
average no, m entrained times maximum pumping rate. The 

number of annually entrained P. hoyi and M. relicta was 

8—1 8 

previously determined to be 1.97 x 10 yr and 1.11 x 10 

yr , respectively, at the Cook Plant (unpublished data, 

GLRD), using a methodology wherein mean densities of 

entrained organisms are calculated for two diel periods (day 

and night) and then multiplied by the proportional amount of 

water pumped by the plant during day and night over the 

sampling interval. Monthly values were summed to estimate 

annual entrainment. For comparison with results of this 

-3 
study, we based the Cook calculation on average number m 

9 3 - 1 
entrained and the maximum pumping rate (3.259 x 10 m yr ) 

at the Cook Plant. These quantities (annual entrainment 

8-1 7 

estimates) were 1.46 x 10 yr for P. hoyi and 8.05 x 10 

yr*" for M. relicta at the Cook Plant. 
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RESULTS AND DISCUSSION 
LAKE MICHIGAN FIELD SURVEY 
General Distributions of Bentho s 

The 10 major constituents of the macrobenthic 
community^ accounting for 109 of 129 taxa (Table 2) and 
99.6% of the average benthic density in the survey area from 
1978-1981 (7,724 m"^) (Table 3), were the Chironomidae (51 
taxa), Naididae (19 taxa), Pisidium (14 taxa), Tubificidae 
(12 taxa). Gastropoda (6 taxa), Turbellaria (4 taxa), 
Enchytraeidae, Pontoporeia hoyi , and Stylodrilus heringianus 
(Table 4). When averaged from 1978-1981, chironomids were 
the most frequently collected of the 10 constituents, 

occurring in 94% of the samples; P. hoyi was the most 

-.9 

numerous (2,726 m ) (Table 3). Other dominant major taxa 

-2 -2 

were chironomids (1,661 m ), naidids (1,256 m ), and 

~2 

tubificids (949 m )• Expressed as a percentage of total 

benthos, distribution of dominant taxa among regions was 
similar when averaged from 1978-1981, with the possible 
exception of P. hoyi (Table 5), Average 4-yr abundance of 
P. hoyi in the inner region was 3,068 m " (42% of mean 
benthic abundance), while that in the outer region was 2,385 
m (29% of mean benthic abundance). 

In 1981, greatest average annual abundance for many 
major taxa was observed in both regions. However, increased 
annual density in the inner region did not occur to the 
extent it did in the outer region. When comparing 1981 
annual abundance (operational year) with that based on 



27 



Table 2. Benthic macroinvertebrates identified from 
Ponar grabs collected from 1977-1981 and from net tows, 
sled tows, entrainment samples, and artificial 
substrates collected during 1981 near the J.H. Campbell 
Plant, eastern Lake Michigan. 

Taxon 



Cnidaria 
Hydrozoa 
Hydroida 
Hydridae 

Hydra sp. 

Platyhelminthes 
Turbellaria (spp. 1-4) 

Annelida 
Oligochaeta 
Prosopora 
Lumbriculidae 

Stylodrilus heringianus 
Plesiopora 
Enchytraeidae (spp. ) 
Naididae 

Amphichaeta leydigii 

Arcteonais lomondi 

Chaetogaster diaphanus 

Chaetogaster diastrophus 

Chaetogaster setosus 

Dero digitata 

Dero pectinata 

Nais barbata 

Nais behningi 

Nais communis 

Nais elinguis 

Nais pardalis 

Nais simplex 

Nais variabilis 

Piguetiella michiganensis 

Pristina f oreli 

Pristina osborni 

Stylaria lacustris 

Uncinais uncinata 

Vejdovskyella intermedia 
Tubif icidae 

Aulodrilus limnobius 

Aulodr ilus pigueti 

Limnodr ilus angustipenis 

Limnodrilus claparedeianus 

Limnodr ilus hoffmeisteri 
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Table 2. Continued. 



Tax on 



Limnodrilus prof undicola 
Limnodrilus spiralis 
Limnodrilus udekemianus 
Peloscolex f rey i 

Peloscolex multisetosus lonqidentus 
Peloscolex superioriensis 
Potamothrix moldaviensis 
Potamothrix ve jdovskyi 
Rhyacodrilus coccineus 
Hirudinea 
Rhynchobdellida 
Glossiphoniidae 

Helobdella stagnalis 

Mollusca 
Gastropoda 
Mesogastropoda 
Hydrobiidae 

Amnicola limosa 

Bithinia tentaculata 

Somatogyrus sp. 
Valvatidae 

Valvata sincera 
Basommatophora 
Lymnaeidae 

Lymnaea sp. 

Physella sp. 
Eulamellibranchia 
Sphaeriidae 

Musculium transversum 

P i s i d i um adamsi 

Pisidium casertanum 

Pisidium compressum 

Pisidium conventus 

Pisidium fallax 

P i s i d i um f erruqineum 

Pisidium henslowanum 

Pisidium idahoense 

Pisidium lill jeborqi 

Pisidium milium 

Pisidium nitidum f . nitidum 

Pisidium nitidum f . pauperculum 

Pisidium subtruncatum 
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Table 2. Continued. 



Taxon 



Pisidium supinum 
Pisidium var labile 
Pisidium walkeri 
Sphaer ium nitidum 
Sphaer ium rhomboideum 
Sphaer ium striatinum 

Arthropoda 
Arachnida 

Acarina (spp.) 
Crustacea 
Amphipoda 
Gamma ri da e 

Gammarus fasciatus 
Gamma r us pseudolimnaeus 
Crangonyx pseudoqracilis 
Haustoriidae 

Pontoporeia hoyi 
Talitridae 

Hyalella azteca 
Mysidacea 
Mysidae 

Mysis relicta 
I sopoda 
Asellidae 

Asellus sp. 
Insecta 
Diptera 
Chironomidae 
Chironominae 
Chironomini 
Chironomus anthracinus -gr ♦ 
Chironomus f luviatilis ~gr . 
Chironomus halophilus -gr > 
Cladopelma sp. 
Cryptochironomus sp. 1 
Cryptochironomus sp. 2 
Cryptochironomus sp.3 
Cryptochironomus cf . rolli 
nr. Cyphomella sp. 
Endochironomus sp. 
Glyptotendipes ( Phytotendipes ) sp, 
nr. Harnischia sp. 
Microtendipes sp. 
Parachironomus cf . abortivus 
Parachironomus sp. 1 
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Table 2. Continued. 



Taxon 



Paracladopelma cf . nereis 

Paracladopelma cf . undine 

Paracladopelma c f . winnelli 

Paratendipes sp. 

Phaenopsectra sp. 

Polypedilum cf. halterale 

Polypedilum cf • illinoense 

Polypedilum cf . scalaenum 

Polypedilum c f . simulan¥ 7 digitif er 

Polypedilum cf • tuberculum 

Polypedilum sp. 2 

Robackia c f . d emei jerei 

Saetheria c f . tylus 
Tanytarsini 

Cladotanytarsus sp. 

Micropsectra sp. 

Paratanytarsus sp. 

Rheotanytarsus sp. 

Tanytarsus sp. 
Orthocladiinae 

Cricotopus (C. ) tremulus ~gr . 

Cricotopus (C.) sp. 

Cricotopus (I^. ) cf. intersectus 

Cricotopus (!_.) cf. suspiciosus ~gr . 

Cricotopus (I_.) sylvestris -gr. 

Cricotopus ( I_. ) sp. 

Heterotrissocladius c f . changi 

Heterotrissocladius c f . bliveri 

Hydrobaenus sp. 

Nanocladius sp. 

Orthocladius (0.) sp. 

Orthocladius (E. ) sp. 

orthocladiini sp.2 

Paracladius sp. 

Parakief f eriell a sp. 

Psectrocladius c f . simulans 

Thienemanniella sp. 
Diamesinae 

Monodiamesa c f , tuberculata 

Potthastia cf . longimanus 
Tanypodinae 

Procladius sp. 

Thienemannymia -gr . 
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Table 2. Continued. 



Taxon 



Ceratopogonidae 
Ceratopogon sp. 
Culicoides sp. 
Probezzia sp. 
Trichoptera 
Molannidae 

Molanna sp. 
Leptoceridae 

Nectopsyche sp. 
Unidentifiable first instar larvae 
Odonata 
Coenagrionidae 
Enallagma sp. 
Hemiptera 
Corixidae 
Pleidae 

Plea striola 
Collembola 



32 



u 

O U 

e C 



I 
<D u 
U >i 
C 

a> o 

V-i CO 
u V- 

U II 

o 
o c 

O '- 

00 

>iCh 

u^ 
c 

0) o 

<u cx> 
u r^ 



c 

CO 



c 

CO 
en 



o o 

U -r-l 



0) 



CM 
I TJ 

e 0) ^ 

4-> CO 

• O »J 
O 0) 
CrH 



o 

>i u 

•r-l W 

en (1) 

CrH 

e -p 

CO C 

en CO 



c 

<u 
-p 

0) 
CO 
0) 



TJ 



c 
(0 
0) C 04 



en rH 
CO a (U 

o a 
u e 
en CO 
o 
o 



o 
a 
a 
< 



CO 



:i: 



6 

o 

Q) C n 
^ O 

CO cOct: 



00 00 

r- 0) 

0) "- 



00 00 



ii>^r--r-if)(DOOin''-0)TfcooJOoor50JCO c 

CDOOOror-o^cococoininr-^iDinoico cd 
CD or-r-tDco-'-coLDin in o) 



Or^CN00I^0)0)O(X)'*'«3-'^^'^0)C7)CS00O O) 

0"^cMCMCDoocoooo)a)0)0)^cor:)t^t^if) oo 



^OOOCO^OCOO^COCDtDt^O^O^CN 00 

0)OOOt^^Lncor)0)«Nif)if)"'-in^in^t^ cn 
cp cor-tDinn^coinin-r- co a 



CDOlDO^^'^^OOrroOCMCM^CN^CO'^lO t^ 
00OOO'^0)'^^C>J'srr^CMCNtf>C4'^(N^O tD 



OOOOOco<Dnh-o)ti)'^0)t^O<D'^oj(r is> 

OOOOincoinncpooin'^roa)inotDc>jO in 
tp ot^t^tD'^'-coinin ^ D cj) 



CD^-'^'^'^CMCDOOOOOO^OCOinoOCOCO '^ 
CMVV (DCJ^inrj-CMinOrOCM CD CM 

CM T- CM -r- i^ 



co-r-T--.;;fcoincot^co(N<i)'r-tpTtc»3CMin-r-in in 

CD ^coo^cMinin-^co -t-ti- t- 

^ r-corj-in-^cM-r-inin rj- t*- 

CO "T- CO •»- ■*- 0) 



lOOOOCOOCOt^ino-^CDCDOinCD^CDlD CO 

in r-'^ocMCMinoocM^-'- t^ tj- 

CM -r- -r- yj 



coO'^O'^coroincDiD^cMcc'^oj^h-^'^ O 

CO V OOCOCOOOOOCOt^ CDV ^ 

CO CMCOr-h-"^ ^J-^ rt 0!) 

CM T- CM ^ I^ 

CMOOC^COOOCDOO'-'^CMtDCD'^CDinrO"^ ^ 

O CM"«^OJC0CDC0t*-l^tD CDV O 

CD CMCM-^CO-^CMCN -r- m 

•^ CM CM -r- ^) 



in 

(0 C3 



~ (0 
E 

E 



fO (0 <D -r- <C 



V) 




3 




C 




(0 








01 




c 






W 


I- 


-D 


0) 





Q) r 


a 


(0 


> 


0) 73 tn 





<0 - 3 


0) 


•0 Q) - 





ax5 

E 


c 
o 



> w 

< 



<0 O c t. 

Q) V -r- +J -o 

<0 — H- > O 

r "D -r- JI r- 

O - JD o >H 

O to 3 C -H 

D)2 K LU LO 



a 

a 

10 (0 



o o 



E 

(0 03 3 
■O T5 - 
O O T3 

a a- 
o > w 

£- O •'- 

•M Q) Q. 

CS 0) 

C3 a 



0) (Q 

a c 

£- 0) i~ a <D c 

ID C (0 10 r- l-i 

JZ - o 



(0 u 
^ 0) 
£- to 



•H T3 <0 (Q 

O 3 £- C 

- > > 

I I. I 



CD L. 

X) (1) 

C SI 

3 -H 

I- O 



(0 
O 

c 

Q) 
15 



33 



Xi rH 




T3 r- 












"D - 








U rH 




C CO 


O CO o -^ o o o 


"^ 






C (0 


O '^ r^ in o ^ CD 


0) 




fl3 -H 


CO T- T- 


in 






(0 "H 


CO -r- 


(P 


fO <D 




£- 










c 






<u x^ 




O -H 










O -H 






- cu 




















C U. rH g 




















••H 0) .rH (d 




















^ -P l-l o 




















-p p a 






















> re 






ro 


CO o O O O O 0) 


.^ 






(0 


CO CM CP ^ O '^ 00 


h- 






4-' 


CNi 


■^ 






■H 


CM ••- 


in 


TJ cr» • 






o 

















TJ c c ^:) 






K 










H- 






0) CO -rH 






















4-> U 0) 






















u V4 p x: 






s 


"^ -^ O O O O r^ 


^ 






00 


t^ CJ) -^ CO O CO !^ 


o 


0) TJ -M 




C 


0) 




'^ 




c 


CJ) 


▼- 


^ 


rH C 




o 


■^ 











'^ 






rH C W J^ 




D) 










D) 








O -H X) (0 




s 










0) 








U CO 0) 




£- 


o 

cx) 


CO IT) O O O O If) 


CO 




J- 


o 

00 


in 0) -r- ^ O ^ ''t 


^ 


<D U C 




L. 


0) 


«r- 


CM 




1- 


0) 


•r- 


CO 


tn ^ C7^ 




0) 










0) 








4J -P C 




3 










4-' 

3 








-rH U (0 




O 










O 








c c ro en 

P -rH C -rH 






0) 
0) 


IP oi O O O O CO 


- 






CJ) 
CJ) 


co t^ O CM o -^ in 


CO 
CM 


O^ 






















U 6 04 u 






















-fH -fH 
























E 




00 






E 




00 






CO 




0) 


O CO O O O O -^ 


'St 


10 




2 


in t^ CO CM o "^ CM 


o 

CO 


C <D 






















O O T> -2«<J 




















X +J <D (0 




















rd -M ^JJ 




















4-> CO u 




















<U C 






















4J 4-) rH U 






^_ 










t— 






C Cd rH 0) 






(0 


CN -^ O "»- O O t»- 

CN 'T- 


^ 






(0 


t^ CM h- CO O CM 00 
CN "r- 


CJ) 

in 


0) O -P 
























>H cn u w 






H 










1- 






(U C (0 






















MH -r-l 0) 0) 






















•rH J^ 0) g 






co 


CO «) O O O O -"a- 


CO 






00 


t^ cp CO CN O CN r- 


r^ 




c 


CJ) 




CM 




c 


CJ) 


f- 


CO 


TJ D ^ O 




o 













"*- 






u u 

>1 u w ^ 




0) 


















rH O 0) 
XJ rH U 






o 

00 
CD 


Tt tp O -^ O O oi 


CO 
CN 






O 

00 
0) 


'r- t^ CO ^ O O "^ 


CP 

CN 


fO (n a (d 




0) 










0) 








-rH ex e <D 




c 
c 










c 








MH P (0 >i 












h-t 








-H O O) 






0) 










CJ) 






-p u ^ 






Cf) 


CM CO O O O O CM 


CM 






0) 


in 0) CO o O O CO 


O 

CO 


c en u 






















0) • (0 






















ro o r- 0) 






















-rH -rH CO 






00 










00 






6 CTi iP 








Tt tp o O O O CM 


CN 

CM 






Cf) 


^ r^ -^ CM o O '^ 


CO 
CN 


MH O '— O 






















o a 1 






















O CO Ui 






















jh X r^ <u 






















0) fd ch ^ 






















X) -P -- O 
















S 4J 




(0 








0) 




P Jh B u 




•D 








X5 

a 




% o o o 




o 
a 










•r-y U 




Q) > 








0) > 




cd HH ro 


c 




(0 0) O 








(0 0) 
•D (0 0) <0 




X 


^ T3 r- T5 








•D »- -D 




^ (n (d 


(0 


E (U ••- E 


Q) O 








E 0) "- E 


0) 




0) ^ O -4J 


H 


O (0 O 3 


a a 








(0 U 3 


a a 






C "D - - 
"- M- XJ 


o 

L. J. C 


^ 






C "D •'- •'- 
- «♦- "D 



L. i. L. 


,-. 


rH -P -rH >i C 




£- XJ - - 


0) -H QJ 


<0 






t. -0 -r- -r- 


^ t: ^ 


(Q 


^ O^rH Cd 




•- — iD W 
H (0 D •■- 











•r- — D W 

r (0 3 "- 







Cd ^ <U P rH 




U Z H Q. 


O C5 O 


H- 






O 2 H CL 


O C3 O 


h- 


b^ O U n Ch 



























34 



C (0 
(0 -H 
L. O 



c 
o 

D) 

L. 

£- 
0) 
+-• 
3 
O 



E 



(U 

a 

-ri 

C 
O 
O 



rH 

(d 



CO T- T- -^00 



CM -r- r^ 



CM in 



coooincoocot^ If) 



OOtDt^Ocoin O 
CM in 



0)Oint»-'«^cMin O) 



CM'^OOO'^CO'SJ- 0) 



^•r-CPCOOCOCM O) 
CM 1- -r- 4n 






OOCOIPCPOCM^ ^ 



"^t^inco-^cMCD O 
CM in 







(n 






T3 













a 


0) 




> 


(0 <D 




o 


•D CO 




0) (Q 


- -D 




r- -0 


E 0) -r- 


E 


(D 


O (0 O 


3 


a a 


C -D - 







o -^ «^ 


•D 


t. t> L. 


I- T) -r- 




Q) 4^ Q) 


•r- -r- n 


m 


n w r 


H 03 3 


•f- 


•H (0 +^ 


O Z H- Q. 


o o o 



C CO 

£- O 

C5 -H 



'«-coo)cMCMcrin oo 

CO -r- "•- -r- 00 





00 


c 


0) 







0) 




0) 




i- 


o 




CO 


;. 


O) 


0) 




c 




c 





0) 






r- 


CO -r- (X) CX) OJ CO T3- 


CM 


0) 


■»- ^ 


in 



OcMooa)"^incM i^ 
CO "- -r- r- 



-r- t^ iS) CD ^ "^ CO in 
rj in 



t^cor^oo-^coin O) 



cnoino^inoD in 



cMr^r^cO'-coin co 
oi in 



rj-T-CD-^CMtDCM in 

CM T- ^ -r- r- 



i^coinr-«0'3-ci^ O 
•r- in 



f^t^inco«^coin ^ 



cotDincDOcoin cd 







■D 

















a 




0) 




> 




ro 0) 




u 




TJ <0 




0) (0 




T3 




r- -D 




E Q) - 


E 


Q) 




(0 


3 


a a 




C "D - 


•r- 







o •- »^ 


•0 


c t. ;_ 


r— 


;» "D - 




4-' 0) 


(Q 


•r- -^ D 


w 


n tn r 


♦J 


r fO 3 




•H (0 -H 





U 2 H- 


CL 


O CJ o 


K 



35 



t- O 



c 
o 

0) 
0) 

£_ 

L, 
0) 

■H 

3 

o 



E 

in 



T3 
<U 
P 
C 

-M 
C 
O 
O 



rH 



CM ^ -r- -r- ^ 00 



rrcMO-^-^inco o 

CM "T- ^ 1- T- t^ 



-r- -r- in 



int^^ro)^coin ^ 



CM -r- (i) 



inO)h>C7)CMCOC0 CO 

•r- in 



cococo<J)cMin'r- 1- 

CM -r- -r- t^ 



inoh-oo'^co'- Tj 
^ 1- in 



tDO)C£)CJ)'^COlO O 

"- in 



r-co(PoocMincD cm 
^ in 



ino)inoo-»-cMtr) (d 



X3 

o 
a 
> 

o 

Q> (0 
Q) O 

a a 
o 

£_£.£. 

s: u) sz 
o o o 



(C (1) 

•^ T3 

E Q) - E 

O R3 O 3 

C U - - 

O - «^ T3 

t. -O - - 

•»- -r- x) tn 

r ro 3 •'- 

o z K a 



o 
o 

r 
a 

(D 



L O 



-^ O) Oi rt ^ iO CO <J) 
in ▼- -r- T- CM CM 



-r-t^CMCM-^inoO <T> 

^ - -- - -^ O 



•^CMCDOCM'^CX) CO 
CO "^ -r- 00 



CJ "^ CD 



rrooo-^coino -^ 

CM '«- T- -r- t^ 



tn'^0)0)cM^r^ !^ 

CM "^ tp 



CMtPOC^CMtPOO CD 

rr ^ ^ ^ ^ O 



C7)^C0C0'^^CD O 
CM -^ -^00 



{£> O t^ G) "T- rt Q) r^ 
CM "^ CD 



cDc^cJ5"^cMinr^ cm 

CM T- -r- t^ 



CDOt^CM-^COt^ CD 
CM •»- -^ CD 







tn 






•D 






O 






a 


0) 




> 


«3 0) 




o 


■D <C 




0) (0 


"- -D 




- TJ 


E Q) -r- 


E 


0) 


O (T3 U 


3 


a a 


C "D "- 







- u- 


-5 


£- L C 


£> -0 - 




<D V 0) 


-r- r- n 


w 


r w i: 


r (0 3 


•r- 


-H (0 -H 


o Z »- 0. 


o o o 



36 



en I 
c o en 

(DUO) 

rH O 

(0 GO OPk 

O CJ> II 



<U 6 
^ O 

4-> V^ 



O no -r-l • 

0) cnc 



<D (0 

u en 



en u 

-P rH---^ 

c »H <u u 



0) o 

o u 
u 

0) (U 



U 0) 



ex U 0) -i^ 



C 
(0 






-^ (0 
CM Q 
I 

e 



u 

(0-—' c 

U 0) 
0) 4J 
4-> en 

P 
• o 

O <UTJ - 

C W C 4-» 

^•H rd 
u 

<D a- 
o c 

U (DfH 
O4-M 0) 

o 0) ex 
>-i u 6 



(d 



c 

(0 



u 
c 

(0 

c 
CO 



fO u 

e e 
o 

r-H C 

o 



CO 

c 
c 



c • 



0) 

CO (UX 
4J^ -P 



in 



u 

CO 

<u 

c 



o g 

(0 S-i 
giM 

0) ^-^ 

X) U CX)I 
(d CO CTl u 



cs 00 
0) -^ 



00 00 
t^ 0) 
CD '- 



cNT-0'^0)OcoO'^'^"^in'^'^-»-'»-oi'»- I 

TtV V-^CO"^"^ V vvvv V 



^v v-^co-^'»-v vvvv V 



in f -^ T- V vvvvvv 



CO -^COCN V VVV 



n n c4 "^ -^ w vvvv 



CO'^OOiOCDOeODCDCDCOCDCJinr^^CN t-- 

ir>v cNOrocNoocpO)Cp id n 

O ^CMDt^ D "^ CO 



0)rooo)0^^tp"^000)t*-cNCNr^r^co a 

rr win'^t^t^ayt^D ip co 

00 ^ m (S> '^ "^ in C4 oj 

CO «r- CM 'T- ^ Q) 



cNOOOoococoin-^^cMOoincocNr-r^ is> 
'^ '*coo)'r-h*if)r^O CO If) 



CO 



O -^ in ^ 



CN 



OOOO'^incooorr^cocooin-^-'-O tc 
D r^o)oo^cMCocoO -^in in 

CN Oi CN r^ 



cNOOOcocDcncNcniDcooirocNOcooo 
in r^ncNf^inoo(X)co o 

CO '^ CO t- Oi -r- 



OC 



>lr3 

riu 

m 

<0 <o 



0) 



0) 

0) "D w 

fl3 - 3 

0) -D <D - 

• (0 (0 <D -^ CO -^ 

D.V -H (0 U £- £- 

tn -r- Q) "D -^ -H T3 

E <0 - «*- > O 

r 13 -- jc — 
c - ' 

O O (0 



15 U . 

h- UJ t/> 



03 
T3 
E 

a 3 

H O — 



u o 






Q) Q) £- _ . _ 

a c (0 w — »-i 
u 



a 
a 

<D U 
^ Q) 



V) 



c 



£-•0(0 
0) 3 £- 



r L. XJ T5 



■H -^ > > 3 -H 
O I X I 






37 



w 

o 

r 

c 



CD 00 
h- 0) 

0) -r- 



00 
0) 



OJ V OiO'^'*- vvvv V 



OO^OOr^CM0)CMTt-^in'r-«r-'r-^C£>^ I 

CO V cmot-^ vvvv V 



O 

00 



CMOOO!^OCMCSr)'!3'"^C0^"^^'«-CM"^ I 

CO oico'^'^ vvvv V 



c 
o 

D) 
(D 

D 

a 



0) 
O) 
(0 

+^ 

c 

0) 

u 

£- 
0) 
CL 



0) 



00 



00 00 
en "^ 



^0'^Or^cor-'^cocM'»-h--«-'r-'»-ooO'^ i 

COV"r-CO"*-"»- VVV V 



OOOOcot^O)in"^cO'^^'^'»-"»-"«-^'»- I 
orf coco-^"^ vvvv V 



inO'^OoJO)coococorrot^inif)0)CMco •»- 

00 Ot^ooh-t^inooo c\ o 

CO CJ)h'"^'-'^CM'»-rJ Tt <r- 

CN ••- Oi "^ •«- 00 



00 
CT) 



coOcoo«>CJ)c>ioo^rj-cMcoc>icoco^int^ 

00 O'^O'^OO'^TJ-CO CM-^ 

O Or^cMCj)'r-TtT-in tp 

CO cs n -r- -^ 



O 



w 

c 

u 

c 

(0 
Q) 



O 
00 
CD 






i^OOOcncpcocD'^coocM'-innt^cDir) o 

to cj^OcMcooocpon-^ -^ rj- 

CM OO'-OOOO'^CNin >*- o 

CM -r- CM t^ 



i^O'^OcoocM-'-Ot^iP^inconO'^'^ rf 

CO co'T-cocMininOiTJ r^v cn 

00 '3-t^coo)CM"»-'»-in w CO 

CN f- CM -^ 00 



CMOOOOOtD'-h-CDO'^Cnt^CDt^d)"- ^ 

If) r^oocMOr^r^ooO oov o 

CO cNincoO'^co CN o 

"^ C\ CM -^ "^ IP 



0) 

c 

•rH 

C 
O 
O 



tn 

0) 

rH 





w 


^ 31 




4-» 


> 


(0 







r 


o 


w 


(0 


(0 


•»- 


M- 



> w 



0) 

E 



in 

C 
CO 

D 

C 

0) 

Q) r 

T5 Q) »- 

Q) -r- TO -^ 

to O £. L. 
Q) "D — -H -D 
(0 - M- > O 



o x: T5 - r r- 

C U •- £1 O > 

O O CO 3 C -H 

i- D5Z K UJ 10 

r — 

a o 



(0 

■0 

o E 
a 3 
o -^ 

y) in 
CO - 
O 0. 



- CO -^ 

Q) 0) L. 

a c CO 

£- X3 CO CO 

O 3 J» IL 



a 
a 

cn CO 

CO o 

• £. CO 

a CO c 

cn .- ►-! 



CO 

o 

c 

0) 
JD 



- > > 3 -M 
O I I X 



0) L 

1} Q} 



38 



1978-1980 (preoperational years), percent increase for 
nearly all 10 major taxa was greater in the outer region 
than in the inner region. The most notable increases were 
those for tubificids, S. herinqianus , and P. hoyi . Although 
all three taxa have displayed consistently greater densities 
in the outer when compared with the inner region during 
previous years, the former two taxa increased 108% and the 
latter 43% during 1981 when compared with 1978-1980. These 
increases were considerably more than concomitant increases 
of approximately 10 to 25% in the inner region when making 
the same comparison. Other representatives of the major 
taxonomic groups displayed rather similar regional changes 
or, as in the case of turbellarians, have exhibited very 
highly variable annual and regional population densities, 
which require interpretive caution. Overall, outer region 
density increased 37% during 1981 when compared with the 
previous 3-yr average, while that in the inner region 
increased only 10%. Nonetheless, density differences 
between 1981 and 1978-1980 were not unusual as annual 
.abundances for all major taxa were quite variable, with 1979 
and 1981 representing years of peak abundance. Subsequent 
evaluation of annual and regional trends will be discussed 
in the analysis of variance (ANOVA) section. 
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Determination of Plant Effect on Benthic 
Populations Due to Discharge of Heated Effluent 

Basis of Interpretation — 

A brief statement is in order to clarify interpretation 
of population changes (R' ) occurring below the limit of 
detection or sensitivity level of the ANOVA (R) . Once 
stated, the R and R' parameters will be presented for 
subsequent taxa with the understanding that an 
interpretation similar to that ensuing below is intended, 
thereby greatly reducing excessive repetition. 

The R value established the minimum level of detection 
at a = 0.05 and P = 0.95 or sensitivity of the ANOVA. When 
estimated R' values were less than expected R values, 
density changes in the operational-regional populations were 
below the minimum detectable limit, thereby suggesting two 
possible interpretations (see METHODS) . First, one could 
assume there was no plant effect. Second, there remains the 
possibility that there was a plant effect, but it was not 
detectable. While it is the second possibility that is of 
most concern, there is no statistical way with present 
datasets to distinguish between no effect and a non- 
detectable effect. 

When the nature of the R values are examined, we feel 
that the R value not only provided a necessary delimiter on 
the ANOVA, but also assured that some quantitative change 
was detectable. By setting a = 0.05 and P = 0.95, we were 
sure of low probabilities for Type I and Type II errors 
(Sokal and Rohlf 1969). By further restricting a and P, we 
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would have decreased the sensitivity of the ANOVA, i.e., 
increased the R value, and have been even more sure that 
excessive R' values would have been indicative of plant 
effect. Alternately, by relaxing a and P to lower levels, 
acceptance of a plant effect would have much greater 
potential for Type I and Type II errors. Consequently, any 
absolute density changes occurring below the limit of 
detection are regarded as having had minimal to no impact on 
the population. In conclusion, we feel R values generated 
at a = 0.05 and P = 0.95 provided sufficiently reasonable 
assurance of detection of a true plant effect while limiting 
Type I and Type II errors, but can not provide absolute 
assurance in a system defined by probability statistics. 
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Pontoporeia hoyi — 

Evaluation of heat effect on P. hoyi density was based 
on the population encountered at 9 to 15 m, as few or only 
sporadic occurrences were observed at depths less than 9 m 
(Fig. 2). Annual monthly abundances followed a similar 
pattern of lowest densities in April and highest in July 
(Fig. 3). A more detailed examination of density trends 
indicated that, although regional density trends were not 
generally widely disparate within a given depth, there was a 
high degree of variability associated with density of 
P. hoyi across depth, month, and year factors (Fig. 4, 
Appendix 1), This variability was expressed by the 
significance of higher-order interaction terms from the 
ANOVA (Table 6) . 

Of main effects in the P. hoyi ANOVA, all except region 
were significant. As was expected, temporal and depth- 
related density differences were evident. Regional density 
differences during each year favored higher numbers of P 
hoyi in the outer when compared with the inner region, 
consequently it was surprising that the regional main effect 
was not significant. Regardless of density differences 
noted by the ANOVA factors, density changes observed for 
P. hoyi at 9 to 15 m were below the level of detection of 
the ANOVA [R' = 1 .22, R = 2.34 (Table 7)], indicating no 
measurable plant effect during the 1978-1981 time period. 
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Table 6. Analysis of variance results for densities 
[log.Q(x+1)] of Pontoporeia hoyi occurring at 9-15 m 
from T978-1981 near the J.H. Campbell Plant, eastern 
Lake Michigan [NS = no significance (p > 0.05), * = 
0.01 < p ^ 0.05, ** = 0.001 < p ^ 0.01, *** = p < 
0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region(R) 


4.10 


1 


4.10 


9.76 


NS 


Depth(D) 


73.34 


2 


36.67 


146.88 


*** 


Month (M) 


86.76 


2 


43.38 


5.38 


* 


Year(y) 


26.63 


3 


8.88 


31.71 


*** 


RD 


0.03 


2 


0.02 


0.03 


NS 


RM 


1.28 


2 


0.64 


0.93 


NS 


DM 


1.87 


4 


0.47 


1.09 


NS 


RY 


1.27 


3 


0.42 


1.50 


NS 


DY 


1.51 


6 


0.25 


0.89 


NS 


MY 


48.36 


6 


8.06 


28.79 


**» 


RDM 


2.00 


4 


0.50 


0.51 


NS 


RDY 


3.65 


6 


0.61 


2.18 


* 


RMY 


4.12 


6 


0.69 


2.46 


* 


DMY 


5.21 


12 


0.43 


1.54 


NS 


RDMY 


11.85 


12 


0.99 


3.54 


*** 


Error 


101.62 


360 


0.28 
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Chironomidae-- 

The average density of chironomids encountered during 
1981 (1,743 m ) was very similar to the average abundance 
of the previous 3 preoperational yr (1,634 m ). Mean 
densities observed at each respective depth (Fig. 5) and 
during each respective month (Fig. 6) during 1981 did not 
differ from the previously observed range of values. Annual 
variability associated with depth and monthly mean densities 
was quite large. Outer region chironomid densities were 
slightly greater than those observed in the inner region at 
most depths, months, and years (Fig. 7, Appendix 2). 
However, based on the chironomid ANOVA, non-significance of 
the regional main effect indicated there was no regional 
density difference for the chironomids. Year, month, and 
depth main effects were significant (Table 8). The very 
high significance of higher-order interactions indicated the 
extreme variability associated with the predictive value of 
main effects variables on the chironomid population. These 
data supported the suspected case among chironomids that 
they are very mobile, transitory animals subject to a 
physically controlled environment in both larval and adult 
stages. Nevertheless, given density fluctuations observed, 
an R' value (1.37) less than the minimum detection limit (R 
= 1.74) (Table 7) indicated there was no detectable plant 
effect on the 3- to 15-m chironomid population during 
1978-1981. 
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Table 8. Analysis of variance results for densities 
[log-,rj(x+1 ) ] of chironomids occurring at 3~15 m from 
1978-1981 near the J.H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0.05), * = 0.01 
< p < 0.05, ** = 0.001 < p ^ 0.01, *** = P ^ 0.001]. 







Sum of 


Deg 


rees of 


Mean 






Parameter 


squares 


f 


reedom 


square 


F-ratio 


Signif . 


Region(R) 


6.11 




1 


6.11 


5.01 


NS 


Depth(D) 


53.66 




4 


13.41 


4.64 


* 


Month(M) 


66.04 




2 


33.02 


6.09 


* 


year(y) 


20.33 




3 


6.78 


33.90 


*** 


RD 


0.92 




4 


0.23 


0.14 


NS 


RM 


1.33 




2 


0.67 


0.66 


NS 


DM 


116.94 




8 


14.62 


5.87 


*** 


RY 


3.65 




3 


1.22 


6. 10 


*** 


DY 


34.67 




12 


2.89 


14.45 


*** 


MY 


32.51 




6 


5.42 


27. 10 


*** 


RDM 


1.93 




8 


0.24 


0.38 


NS 


RDY 


19.17 




12 


1.60 


8.00 


*** 



RMY 6.08 6 1.01 5.05 *** 

DMY 59.70 24 2.49 12.45 *** 

RDMY 15.10 24 0.63 3.15 *** 

Error 119.46 600 0.20 
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Naididae — 

Annual naidid abundance with respect to depth was 
highly variable. Greatest densities were most consistently 
encountered at 9 to 12 m, but occasional large abundances at 
6 and 15 m contributed much to observed variability. Only 
at 3 m were consistently similar annual densities noted 
(Fig. 8). Monthly naidid abundances during each year tended 
to be greatest in July, but were quite variable across years 
(Fig. 9). Despite these differences, 1981 mean densities by 
month and depth were not dissimilar to the range of values 
observed previously. While mean naidid abundance increased 
17% during 1981 (1,408 m ) when compared with a similar 
estimate from 1978-1981 (1,206 m ), proportionate increases 
in each region were similar. Regional comparisons at each 
depth indicated greatest discrepancies occurred at 3 and 6 
m, with remaining depths exhibiting very similar abundances 
and trends. Regional differences at 3 and 6 m were noted 
entirely in July. Annual regional differences at 6 m were 
sporadic and lacked consistency. Similar estimates at 3 m 
indicated the occurrence of a consistently occurring 
population in the inner region,- but fluctuating numbers in 
the outer region (Fig. 10, Appendix 3). 

The naidid ANOVA indicated very highly significant 
year, month, and depth main effects as well as interaction 
terms. However, the region main effect was not significant 
(Table 9). The R' value (1.30) for the naidid population 
occurring at 3 to 15 m during 1978-1981 was considerably 
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Table 9. Analysis of variance results for densities 
[log-Q(x+1)] of naidids occurring at 3-15 m from 
1978-T981 near the J.H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0*05), * = 0.01 
< p ^ 0.05, ** = 0.001 < p < 0.01, *** = P ^ 0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region (R) 


2.64 


1 


2.64 


17.60 


** 


Depth (D) 


410.22 


4 


102.55 


24.30 


*** 


Month{M) 


436.57 


2 


218.28 


62.54 


*** 


YearCY) 


19.22 


3 


6.41 


13.64 


*** 


RD 


13.52 


4 


3.38 


2.82 


NS 


RM 


23.25 


2 


11.62 


8.24 


* 


DM 


70.61 


8 


8.83 


3.25 


* 


RY 


0.45 


3 


0.15 


0.32 


NS 


DY 


50.60 


12 


4.22 


8.98 


*** 


MY 


20.95 


6 


3.49 


7.43 


*** 


RDM 


20.94 


8 


2.62 


2.32 


NS 


RDY 


14.34 


12 


1.20 


2.55 


** 


RMY 


8.46 


6 


1.41 


3.00 


** 


DMY 


65.36 


24 


2.72 


5.79 


*** 


RDMY 


27.17 


24 


1.13 


2.40 


*** 


Error 


280.22 


600 


0.47 
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lower than the R value (2.33) (Table 7), indicating no 
measurable plant effect was associated with naidid density 
fluctuations. 
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Tubif icidae-- 

Mean abundance of tubificids during 1981 increased 108% 

-2 
from the 1978-1980 preoperational average density of 749 m 

-2 

to 1,547 m . Annual tubificid densities during 

preoperational years were quite similar with a range of only 

-2 
627 to 886 m . Although tubificid populations tended to be 

quite patchy and contagiously distributed, the homogeneous, 

fine sandy substrate type which generally lacked significant 

quantities of silt may have influenced the observed low 

annual densities and variability of tubificids. 

Nonetheless, within certain depths and months annual 

variability was quite high, particularly at 12 and 15 m 

(Fig. 11), and during July (Fig. 12). Results from the 

tubificid ANOVA based on population densities occurring at 9 

to 15 m indicated significant annual and monthly density 

differences, no significant difference among depths, but 

highly significant higher-order interactions (Table 10). 

Generally, within the 9- to 15-m depths, abundance of 

tubificids was quite similar (Fig. 11). A significant 

difference among regional tubificid abundances confirmed the 

previously observed trend of outer region tubificid 

densities being greater than those in the inner region. 

This trend was observed consistently at all depths. While 

the absolute number of tubificids in each region differed, 

density trends were generally similar even though the trend 

in the outer region was more pronounced, particularly at 9 

and 12 m (Fig. 13, Appendix 3)o No definitive explanation 
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Table 10. Analysis of variance results for densities 
[log-Q(x+1)] of tubificids occurring at 9-15 m from 
1978-T981 near the J.H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0.05), * = 0.01 
< p < 0.05, ** = 0.001 < p ^ 0.01, *** = P ^ 0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region(R) 


18.06 


1 


18.06 


13.18 


* 


Depth(D) 


7.19 


2 


3.60 


1.46 


NS 


Month(M) 


19.55 


2 


9.77 


5.65 


* 


Year(Y) 


20.93 


3 


6.98 


15.17 


*** 


RD 


2.53 


2 


1.27 


0.76 


NS 


RM 


0.41 


2 


0.20 


0.50 


NS 


DM 


3.09 


4 


0.77 


1.24 


NS 


RY 


4.10 


3 


1.37 


2.98 


* 


DY 


14.79 


6 


2.46 


5.35 


*** 


MY 


10.40 


6 


1.73 


3.76 


*** 


RDM 


1.27 


4 


0.32 


0.23 


NS 


RDY 


10. 10 


6 


1.68 


3.65 


** 


RMY 


2.41 


6 


0.40 


0.87 


NS 


DMY 


7.42 


12 


0.62 


1.35 


NS 


RDMY 


16.86 


12 


1.40 


3.04 


*** 


Error 


167.21 


360 


0.46 
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can be offered to explain the regional difference, except 
that the difference was a regional dissimilarity inherent at 
the start of the study. Nonetheless, it was evident that 
despite this initial regional difference, no measurable heat 
effect due to plant operation was associated with density 
changes observed in the tubificid population at 9-15 m 
during 1978-1981, since the R' value (K54) was much lower 
than the R value (2.98) (Table 7). 
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Enchytraeidae — 

The 1981 average number of enchytraeids (128 m ) was 
the same as the 1978-1980 preoperational mean density. 
Likewise, mean density at each depth and average monthly 
abundance during 1981 were quite similar to the previously 
observed range of values. Greatest enchytraeid abundance 
was observed at 12 and 15 m (Fig. 14). As very few 
enchytraeids were encountered at depths less than 9 m, the 
enchytraeid ANOVA considered only the 9- to 15-m population 
densities. Although there was little difference among 
monthly density estimates for enchytraeids, there was a 
slight increase from April to July and October, which had 
fairly similar enchytraeid densities (Fig. 15). 

Based on the enchytraeid ANOVA, there were significant 
year, depth, and regional main effects, no significant 
monthly differences, and variable significance among higher- 
order interactions for population densities occurring at 9 
to 15 m (Table 11). Examination of regional enchytraeid 
density trends at 9 to 15 m suggested only occasional strong 
regional differences, primarily during October at 9 and 12 m 
and April at 15 m (Fig. 16, Appendix 3). However, overall 
operational-regional density differences were negligible. 
As the R' value (2.81) was below the detection limit of the 
ANOVA (R = 3.81) (Table 7), no measurable plant effect was 
associated with the enchytraeid population occurring at 9 to 
15 m from 1978-1981 near the Ccimpbeli Plant. 
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Table 11, Analysis of variance results for densities 
[log-Q(x+1)] of enchytraeids occurring at 9-15 m from 
1978-1981 near the J.H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0.05), * = 0.01 
< p < 0.05, ** = 0.001 < p < 0.01, *** = P ^ 0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region (R) 


27.48 


1 


27.48 


12.55 


* 


Depth (D) 


74.90 


2 


37.48 


10.12 


* 


Month(M) 


48.47 


2 


24.23 


2.64 


NS 


Year(Y) 


71.17 


3 


23.72 


33.89 


*** 


RD 


1.96 


2 


0.98 


0.80 


NS 


RM 


2.56 


2 


1.28 


0.56 


NS 


DM 


14.26 


4 


3.57 


2.63 


NS 


RY 


6.57 


3 


2.19 


3.13 


* 


DY 


22.20 


6 


3.70 


5.29 


*** 


MY 


55.07 


6 


9.18 


13.11 


*** 


RDM 


13.45 


4 


3.36 


1.80 


NS 


RDY 


7.37 


6 


1.23 


1.76 


NS 


RMY 


13.71 


6 


2.29 


3.27 


** 


DMY 


16.35 


12 


1.36 


1.94 


* 


RDMY 


22.42 


12 


1.87 


2.67 


** 


Error 


252.14 


360 


0.70 
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Stylodrilus heringianus — 

During 1981, as was observed during the previous years, 
S. heringianus was found consistently only at 15 m 
(Fig. 17). A large increase in the mean density of 
S. heringianus during 1981 (252 m ) compared with the 
1978-1980 preoperational average density (128 m ) resulted 
from an unusually high number of individuals encountered 
during July, 1981 (Fig. 18). When comparing 1981 regional 
abundances of S. heringianus with the 1978-1980 regional 
averages, mean density increased 108% in the outer region, 
but only 14% in the inner region. This regional difference 
was mostly attributable to densities encountered in the 
outer region at 15 m during July (Fig. 19, Appendix 3). In 
addition, there continued to be a regional density 
difference in October at the same depth. While this trend 
clearly supported an apparent regional difference favoring 
higher abundances of S. heringianus in the outer region than 
in the inner region, the regional main effect in the ANOVA 
was not significant (Table 12). Both year and month effects 
were significant. The value of R (13.35) (Table 7) was very 
high by comparison with those determined for other taxa, but 
likely reflected the high degree of variability one would 
expect of a population occurring at the nearshore extent of 
the main body of the population. A similar calculation of R 
at the Cook Plant for S. heringianus occurring in an 
increasingly more favorable habitat at 16 to 24 m was 6.70 
(unpublished data, GLRD). The R' value of 1.24 indicated 
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Table 12. Analysis of variance results for densities 
[Ioq-qCx+I)] of Stylodrilus herinqianus occurring at 
15 m from 1978-1981 near the J.H. Campbell Plant, 
eastern Lake Michigan [NS = no significance (p > 
0.05), * = 0.01 < p ^ 0.05, ** = 0.001 < p ^ 0.01, 
*** = p ^ 0.001]. 







Sum of 


Deg 


rees of 


Mean 






Parameter 


squares 


f 


reedom 


square 


F-ratio 


Signif . 


Region{R) 


12.17 




1 


12.17 


3.58 


NS 


Month(M) 


27.08 




2 


13.54 


13.41 


* 


year(Y) 


23.99 




3 


8.00 


9.09 


*** 


RM 


1.08 




2 


0.54 


0.24 


NS 


RY 


10.20 




3 


3.40 


3.86 


** 


MY 


6.07 




6 


1.01 


1.15 


NS 


RMY 


13.56 




6 


2.26 


2.57 


* 


Error 


105.24 




120 


0.88 
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that from 1978 to 1981 there was no measurable plant effect 
associated with the S, herinqianus population occurring at 
15 m near the Campbell Plant. 
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Gastropoda — 

The mean number of gastropods encountered during 1981 
(156 m ) was considerably greater than a similar estimate 
from the previous 3 yr (56 m )• However, regional 
gastropod increases were fairly similar, with a 33% increase 
in the inner region and a 54% increase in the outer region. 
Very few gastropods were collected at depths less than 9m. 
Greatest gastropod densities were observed at 12 and 15 m 
(Fig. 20). Monthly gastropod densities tended to increase 
from April through October, but were quite variable from 
year-to-year (Fig. 21). During 1981, the higher gastropod 
abundances observed at 12 and 15 m and during July and 
October did not greatly exceed the previous range of values. 
As noted during preoperational years, and as was the case 
during 1981, there were very few regional density 
differences among the gastropods encountered at 9 to 15 m 
(Fig. 22, Appendix 4). 

Based on the gastropod ANOVA, regional and month main 
effects were not significant, while main effects of depth 
and year were significant (Table 13). The ANOVA had an R 
value of 4.46 (Table 7). As the R' value was 1.78, no 
detectable plant effect was associated with gastropods 
occurring at 9 to 15m from 1978 to 1981 near the Campbell 
Plant. 



90 



lO 



o 

O 
Cl 

O 

i_ 

(/) 

o 
CD 



00 (Ti O — 
N- N- OO 00 
(J5 (T> (J) OD 



I I 



O 

o 




( lAJ ON) AllSN^a NV3^^ 



91 






D 
■D 
O 
Q. 

O 

^_ 

C/) 
D 



o> o — 

h. OO 00 
0) 0) (7> 




r 
o 
o 



-T- 
o 
o 

CM 



1^ 

o 
o 



{ ui ou) AijsusQ uoaw 





<U Xi 


o 

00 


^ u 


91 » 


(0 (0 


•§ 


^J J-i 0) 




o 


fe ° 


C M-l C 




TJ V-i -H 




0) (u in^ 


00 


4-J 4-> <D -P 


2> 


U to 4-3 -iH 




<D (0 (0 5: • 




rH (u e J^ 




»H •r^ cn (0 




O C 4->X X> 




O -r-l W 4J 




0) UrH 




in r- 0) (0 




TJ CO >iTJ U 


_ 


O 0^ 4-> •^^ 


00 

22 


0,r- .M rH 4-> 




O in rH u 


o 


U X C (0 0) 


00 


^ 4J en 0) > 


0} 


^ in p p u 


"5 


g fO O 0) >i 


f "^ 


1 CT» J-4 > XI 


2 


X • o 




IM 4-> 4-) TD 


00 


O C CP CD 


fe 


CX) (0 C -P 




--- r- rH -r-l O 




(N ch 04 en c 




I '^ (0 0) 




g rH U TD 




U rH 0) 




V-i 0) <D > >-< 




0) X5 J5 CO O 




Xi O CX u 




6 4J e >i v^ 




D U (0 X> 0) 




coo 




— ' T> TD 




n:) • 0) ui 




>i c :r: 4J (d 




-P ro D T3 


00 
(J) 


.-1 • ac 




in -J^ E m 




C >i O -P 


2 


0) rH 0) U to 


2 _ 


TJ P x: 


•;3 


«-D -P 0) 


?i < 


C V4 • 


g> 


fO ^ >-i (U-^ 




0) r-l CO ^ O 


00 


:s .n4 0) vD 


N- 


ui c 4:: 


O) 


a -p II 




• < c c 




T- CO C 




cNj cp en g -^ 




c ••-» 




. .rH ^ 4:: U 




en V4 u u CO 




.rH P.rH CO <U 




1x4 T5 S 0) >^ 



92 




E 



a 
o 

Q. 

o 

^^ 
-•— 

O 

C9 




± o 
I 




o 



o 



(^-uj 'ou) XiisufliQ uoay^ 



00 
0) 


C 




tn >i 0) 




T5 X> m 




O 0) 


o 


a (U T> u 


00 


O ^ 0) Ou 




u ro -P 


JO 


4-i ►J O U 


o 


tn c (0 


o 

o> o 


m C 0) 0) 


a^ U ro c 


0> 


0) 




M-i 4-> U (0 




o en o 0) 




CO J-4 J-i • 


00 


--- 0) Ui CO fO 


0> 


CM <D <U 


~. 


1 g 4J U 




e O T3 C CO 




V4 }^ 0) 




V4 ^ CO G <D 




0) TJ -P U 




.ri r- c CO c 




£ cr> CO 0) (U 




p cr> 4-> V4 u 




C 'r- C/) 4-> <D 


CO 


— ' u-i 


0) 


^ O 0) 




in CT> • 4-j u 




<U D -U 




•r-i o c en o 


o 

00 


4-> U CO T5 4-) 


C7> 


•H ^ rH C 




w -p cu o en 


— 


c aTJ 


-3 


(L> 00 rH U) C 


(Ti 


T) r- r-i <u o 




m 0) J^ ex 




c ^ ^ j^ en 




CO a o 0) 




<u >^ g O J^ 


00 


S 0) CO V4 


N- 


^ CJ c o 


0) 


rH O O U 




CO -P • -rH 




c u p: en c 




O O (DO 




•rH • V^ -rH 




cnno »^ D^ 




(DC U (U 




Ui CO <D d) V4 




eC C 


_ 


V-i ^ 4J c u 


CO 
0) 


0) >1 hH 0) 




4-) i-H U 4J 




D P CO D 




o n <u • o 


o 


c---^ 


00 


no ^ vD 


(Ti 


c rH e 


— 


CO •'H II 0) 


Ol 


j^ in en 


< 


u a^ c >^ 




0) < 1 "-^ CO 


O) 


a o\ x: 




a c u o 




^-H -rS 4-> CO tn 




CO ^ -H 


00 


T) TJ 




• CD C rH 




CN 4J CO CO rH 




csj u en o CO 




<U .r-l .rH g 




• rH ^ 4-) U 




CT>rH O U <U 




.rH O -H X: 




Cl, u S > -P 



93 




o 

-00 






E 

CM 



D 
TD 
O 
Cl 

O 

^_ 

(J) 
O 





o 
o 

CSJ 



1^ 

o 
o 



o 
o 

ro 



(-UJ ou) AiisudQ uoa^ 



(U 

D 
C 



C 

o 
o 



04 



cn 

-r-» 
fX4 



94 



t-sA 




o 
_oo 



0^ 



00 



E 
in 



D 
■D 
O 
Q. 
O 

CO 
D 




o 




o 
o 



o 
o 



o 
o 



^5 
0) 

"^ C 

< .H 

4-> 
C 

o 
o 



CM 
(NJ 



(^.UJ 'OU) A4JSU8G UD8^ 



•rH 
tX4 



95 



Table 13. Analysis of variance results for densities 
[Ioq-qCx+I)] of gastropods occurring at 9-15 m from 
1978-1981 near the J.H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0.05), * = 0.01 
< p < 0.05, ** = 0.001 < p ^ 0.01, *** = P ^ 0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region (R) 


3.92 


1 


3.92 


3.84 


NS 


Depth (D) 


103.99 


2 


51.99 


76.46 


*** 


Month (M) 


43.63 


2 


21.82 


2.26 


NS 


year(Y) 


18.56 


3 


6.19 


8.72 


*** 


RD 


0.66 


2 


0.33 


0.18 


NS 


RM 


4.03 


2 


2.01 


1.34 


NS 


DM 


6.56 


4 


1.64 


0.88 


NS 


RY 


3.07 


3 


1.02 


1.44 


NS 


DY 


4.09 


6 


0.68 


0.96 


NS 


MY 


57.87 


6 


9.64 


13.58 


*** 


RDM 


0.10 


4 


0.02 


0.04 


NS 


RDY 


11.24 


6 


1.87 


2.63 


* 


RMY 


9.02 


6 


1.50 


2.11 


* 


DMY 


22.33 


12 


1.86 


2.62 


** 


RDMY 


6.96 


12 


0.58 


0.82 


NS 


Error 


254.34 


360 


0.71 







96 



Pisidium — 

Although average Pisidium abundance was slightly higher 
during 1981 (536 m"^) than during the 1978 to 1980 period 
(387 m ), proportionate increases were similar m both the 
inner (17%) and outer (15%) regions. Based on the range of 
1978 to 1980 densities, pisidia abundances observed during 
1981 at each depth (Fig. 23) and during each month (Fig. 24) 
were very similar within respective depths and months. In 
addition, with the exception of 9-m July pisidia densities 
within which outer region densities exceeded inner region 
densities, no consistent regional density trends were 
evident (Fig. 25, Appendix 4). 

The Pisidium ANOVA was based on the 9- to 15-m 
population densities, since few pisidia occurred at depths 
shallower than 9 m. The main effect due to depth was highly 
significant as pisidia abundance steadily increased with 
depth (Table 14). While the year main effect was also 
highly significant, neither month nor regional main effects 
were significant. Across months there tended to be a slight 
increase in average density from April to October, but 
annual variability associated with monthly densities was 
quite high, thereby obscuring a definitive trend. As the R 
value (3.07) was considerably greater than the R' value 
(1.87) (Table 7), we concluded there was no detectable plant 
effect on the pisidia population at 9 to 15 m during 1978 to 
1981 near the Campbell Plant. 
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Table 14. Analysis of variance results for densities 
[log-Q(x+1)] of Pisidium occurring at 9-15 m from 
1978-T981 near the J*H. Campbell Plant, eastern Lake 
Michigan [NS = no significance (p > 0.05), * = 0.01 
< p < 0.05, ** = 0.001 < p ^ 0.01, *** = P ^ 0.001]. 



Sum of Degrees of Mean 
Parameter squares freedom square F-ratio Signif. 



Region(R) 


5.11 


Depth{D) 


158.98 


Month{M) 


17.79 


Year(y) 


7.91 


RD 


3.58 


RM 


4.10 


DM 


2.30 


RY 


4.45 


DY 


30.30 


MY 


15.44 


RDM 


6.53 


RDY 


4.91 


RMY 


3.89 


DMY 


11.19 


RDMY 


8.88 



1 


5.11 


3.45 


NS 


2 


79.49 


15.74 


** 


2 


8.90 


3.46 


NS 


3 


2.64 


5.39 


** 


2 


1.79 


2.18 


NS 


2 


2.05 


3.15 


NS 


4 


0.58 


0.62 


NS 


3 


1.48 


3.02 


* 


6 


5.05 


10.31 


*** 


6 


2.57 


5.24 


*** 


4 


1.63 


2.20 


NS 


6 


0.82 


1.67 


NS 


6 


0.65 


1.33 


NS 


12 


0.93 


1.90 


* 


12 


0.74 


1.51 


NS 



Error 177.64 360 0.49 
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Turbellaria — 

Of the 10 taxonomic groups encountered, turbellarian 
densities were the most unpredictable from year-to-year. 
Neither average densities within depths (Fig. 26) nor months 
(Fig. 27) were predictable on an annual basis. During 1981, 
the proportionate increase in turbellarian density, when 
compared with that of the previous 3 yr, was considerably 
higher in the outer (72%) than in the inner (27%) region. 
Overall, turbellarian density increased 81% during 1981 when 
compared with the 3-yr preoperational average. Examination 
of regional density trends at each depth indicated 
considerable annual variability at all depths, but only 
sporadic regional density differences at the 6- to 15-m 
depths (Fig. 28, Appendix 1). At 3 m, outer region 
turbellarian densities generally exceeded inner region 
abundances, but nonetheless mean densities in both regions 
were highly variable and unpredictable from year-to-year. 

All main effects except month were significant based on 
the ANOVA of turbellarian densities occurring at 3 to 15 m 
(Table 15). Year was the most highly significant of the 
main effects, followed by depth and region. Significance of 
the region main effect indicated the outer region density 
was greater than the inner region density. This regional 
density difference was likely attributable to the large and 
highly variable turbellarian abundance associated with the 
3-m depth regime. With the ANOVA' s sensitivity determined 
to be R = 2.49 (Table 7), disparate regional density trends 
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Table 15. Analysis of variance results for densities 
[log-Q(x+1)3 of turbellarians occurring at 3-15 m 
from 1978-1981 near the J.H. Campbell Plant, eastern 
Lake Michigan [NS = no significance (p > 0.05), * = 
0.01 < p ^ 0.05, ** = 0.001 < p ^ 0.01, *** = p ^ 
0.001]. 

Sum of Degrees of Mean 
Parameter squares freedom square F-ratio Signif. 



Region(R) 22.95 1 22.95 25.50 * 



** 



Depth(D) 69.66 4 17.42 7.00 

Month(M) 144.00 2 72.00 4.12 NS 

year(y) 189.54 3 63.18 117.00 *** 

RD 53.35 4 13.34 4.78 * 

RM 9.85 2 4.92 4.73 NS 

DM 28.26 8 3.53 1.35 NS 

RY 2.70 3 0.90 1.67 NS 

DY 29.88 12 2.49 4.61 *** 

MY 104.75 6 17.46 32.33 *** 

RDM 19.51 8 2.44 1.10 NS 

RDY 33.49 12 2.79 5.17 *** 

RMY 6.21 6 1.04 1.93 NS 

DMY 63.00 24 2.62 4.85 *** 

RDMY 53.11 24 2.21 4.09 *** 

Error 325.97 600 0.54 
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did not result in cumulative differences among regions 
during subsequent years. When contrasting before and after 
operation regional averaged log densities, an R' value of 
U48 indicated no measurable heat effect due to plant 
operation was attributable to density changes observed among 
the turbellarian populations occurring at 3 to 15 m during 
1978 to 1981 near the Campbell Plant. 
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Total Benthos — 

Average benthic density in 1981 (9,715 m ) increased 
38% from the 1978-1980 mean density (7,061 m""^). Although 
1981 benthic densities were generally higher than similar 
estimates from the previous 3 yr, there was little 
difference among preoperational average benthic densities 
observed within each respective depth (Fig. 29) and during 
each respective month (Fig. 30). During 1981, outer region 
mean benthic abundance increased 37% compared with only 10% 
in the inner region, primarily due to increased oligochaete 
and Pontoporeia hoyi abundance in the outer region. 
Examination of regional density trends indicated only 
occasional density differences among regions (Fig. 31, 
Appendix 1), none of which were sustained beyond the 
particular month within which they were observed. 

The ANOVA based on total benthic densities at 3 to 15 m 
indicated all main effects were highly to very highly 
significant except a non-significant regional main effect 
(Table 16). Nearly all higher-order interactions were very 
highly significant. Subsequent contrasting of before and 
after operation regional averaged log densities indicated 
actual population changes were equivalent to R' = 1.37 
(Table 7). As the sensitivity of the ANOVA (R = 1.60) was 
greater than density changes observed, no measurable plant 
effect was evident for the benthic population occurring at 3 
to 15 m in the vicinity of the Campbell Plant from 1978 to 
1981. 
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Table 16. Analysis of variance results for total 
benthic densities [log^fxCx+l ) ] occurring at 3-15 m 
from 1978-1981 near th4 J.H. Campbell Plant, eastern 
Lake Michigan [NS = no significance (p > 0.05), * = 
0.01 < p ^ 0.05, ** = 0.001 < p ^ 0.01, *** = p ^ 
0.001]. 







Sum of 


Degrees of 


Mean 






Parameter 


squares 


freedom 


square 


F-ratio 


Signif . 


Region (R) 


4.63 


1 


4.63 


3.62 


NS 


Depth(D) 


212.39 


4 


53.10 


42.14 


*** 


Month(M) 


128.90 


2 


64.45 


15.53 


** 


Year(Y) 


17.30 


3 


5.77 


41.21 


*** 


RD 


7.01 


4 


1.75 


0.98 


NS 


RM 


4.09 


2 


2.04 


2.62 


NS 


DM 


48.93 


8 


6.12 


4.86 


** 


RY 


3.83 


3 


1.28 


9.14 


*** 


DY 


15.10 


12 


1.26 


9.00 


*** 


MY 


24.88 


6 


4.15 


29.64 


*** 


RDM 


3.16 


8 


0.39 


0.48 


NS 


RDY 


21.49 


12 


1.79 


12.79 


*** 


RMY 


4.65 


6 


0.78 


5.57 


*** 


DMY 


30.19 


24 


1.26 


9.00 


*** 


RDMY 


19.48 


24 


0.81 


5.79 


*** 


Error 


86.75 


600 


0.14 
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Conclusions-- 

Based on the 4-yr design and R' and R values, we 
concluded that no measurable plant effect due to discharge 
of heated effluent into Lake Michigan was observed in 
benthic density fluctuations among nine taxa and total 
benthos occurring in the 3- to 15~m depth range during 
1978-1981. As we have indicated, the nature of the test did 
not totally eliminate the possibility there may have been a 
plant effect below the limit of detection. However, density 
differences below the limit of detection were considered 
negligible, and given the variability in the system, would 
be very difficult to document without greatly increasing the 
sampling effort. Taking into account inherent variability 
and regional differences, density changes exceeding the 
limit of detection would have strongly indicated a plant 
effect, although not absolutely, as there may have been 
other unmeasured factors causing a difference. 

Values of R' were generally quite low, ranging from 
1.22 to 2.81. Low R' values indicated there was overall 
minimal density change among respective population 
densities, in spite of large, sporadic temporal density 
differences. While temporal density differences may be of 
importance for studies of life history, in an impact study, 
these differences are of interest, but by necessity, must be 
viewed on a long-term basis to assess their importance. For 
this reason, it is important to reiterate the high priority 
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to which impact studies must place on long-term data 
collection procedures employing extensive sampling effort in 
the context of a well-conceived statistical design. 
Unfortunately, in the present case, the vagaries of 
constructing and initiating operation of Unit 3 did not 
correspond as planned with the statistical design, resulting 
in an unbalanced design in need of alteration. Clearly, 
more than 1 operational yr would have been preferable. 
While we can only speculate on what effect additional 
operational years may have on our conclusions, we feel the 
alteration of the original design did not have a negative 
effect on results obtained or conclusions drawn. Based on 
our experience at the Cook Plant where results of both 
designs were compared (unpublished data, GLRD) , conclusions 
were identical. 



ARTIFICIAL SUBSTRATE STUDY ON THE INTAKE RIPRAP 

Benthic macroinvertebrates collected from the five 
artificial substrate replicates were represented by 15 taxa 
(Table 17). Chironomids were most diverse being represented 
by seven taxa. Of the seven chironomid taxa, as well as 
with most remaining taxa collected from the artificial 
substrates, all except Paratanytarsus sp. and Rheotanytarsus 
sp. were collected by Ponar grabs from surrounding sandy 
substrates in the lake proper. As both Paratanytarsus 
sp. and Rheotanytarsus sp. form tubes adhering to firm 
substrates, e.g., rock and plant stems, their presence on 
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the riprap was not unexpected. In addition, the abundance 
of Gamma r us spp. ( Gamma r us fasciatus and G. pseudolimnaeus ) 
and Asellus sp. was considerably different than that 
encountered from surrounding sandy substrates. Whereas 
Gamma r us spp. averaged 139 m (8% of total benthic density) 
and Asellus sp. averaged 313 m (19% of total benthic 
density) on the riprap, in sandy substrates neither animal 
was collected prior to 1981. During 1981, occurrence of 
both animals was restricted to the inner region [ Gamma r us 
fasciatus (2.7 m ), Asellus sp. (8.8 m )]. Considering 
the difference between sampling techniques, a fairly similar 
ratio of Asellus -to- Gammarus was observed (riprap substrate 
= 3.3, inner region sandy substrate = 2.2). 

Hydra sp. had the highest density among macro- 
invertebrates occurring on the riprap. Colonies of Hydra 
sp. are suspected to occur over large areas of the riprap 
where water flow and food supply are sufficient. These 
animals occur in colonies of budding, but not separate 
individuals. When samples are processed for sorting 
purposes, there is no accurate way to determine the number 
of individuals present as the washing and handling process 
breaks juvenile individuals prematurely from the adult, 
making it virtually impossible to accurately discriminate 
free-living adult individuals from developing progeny. 

Of enumerated taxa, the most abundant form was naidids 
(38% of the benthos) which was comprised almost entirely by 
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Chaetoqaster diaphanus (99.4%). The only other naidid 
encountered was Nais variabilis . 

The second-most abundant benthic invertebrate on the 
artificial substrates was the turbellarians (26% of 
benthos). Like their counterpart in the benthic lake 
survey, turbellarian densities were quite variable among 
replicates. Other taxa of importance, expressed as a 
percentage of total benthic density, were Asellus sp. (19%), 
Gamma r us sp. (8%), and the snail, Physella sp. (7%). 
Chironomids contributed very little to total benthic 
abundance (1.3%). While it was surprising that the 
chironomids were not more numerous, low densities may have 
been a function of poor sieve retention of early instars 
characteristic of fall months, presence of large numbers of 
Hydra sp., physical location of substrates, and substrates 
themselves. Clearly, for chironomids and other taxa, 
densities would be better estimated by more frequent 
estimates equally dispersed throughout the year. At 
present, data for both density and diversity largely 
represent an autumn estimate. 

The density of benthos on the riprap is important in 
connection with potentially entrainable invertebrates. The 
most likely macroinvertebrates to be entrained are those 
that exhibit pelagic, migratory behavior. Of taxa 
encountered on artificial substrates, amphipods, 
chironomids, and naidids demonstrate the highest degree of 
migratory activity (Marzolf 1965; Wells 1968; Wiley and 
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Mozley 1978). Consequently, these same animals would be 
expected to occur in the water column of Lake Michigan and 
subsequently in entrainment samples. 



EFFECTS OF ENTRAINMENT ON MALACOSTRACAN 
POPULATION BIOLOGY AND LAKE ECOLOGY 

Malacostracan Abundance in the Vertical 
Water Column of Lake Michigan 

General Distribution — 

The percent compositions of malacostracans collected in 

both net and sled tows from Lake Michigan were similar. In 

order of decreasing density, Pontoporeia hoyi , Mysis 

relicta , and Gamma r us spp. (G. fasciatus and 

G. pseudolimnaeus ) were the dominant taxa in samples 

collected by both methods, with the former comprising 93% of 

the malacostracans in net tows and 75% in sled tows. M. 

relicta and Gamma r us spp. comprised a higher portion of 

total malacostracan density in sled tows {14% and 11%, 

respectively) than among net tows (5% and 2%, respectively). 

All remaining malacostracan species, Hyalella azteca , 

Crangonyx pseudogracilis , and Asellus sp. , were observed in 

relative abundances of less than 1% in net and sled tows 

(Table 18). Of these rare malacostracans, Asellus sp. did 

not occur in net tows as it is not a pelagic form, but a 

cryptic omnivore foraging on food sources supplied by the 

riprap substrate. 
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-3 ""3 
Table 18, Average density (no. x 10 m ) for 

malacostracans collected in net tows, sled tows, and 

entrainment samples (Entr) during 1981 near the 

J.H. Campbell Plant, eastern Lake Michigan (n = number of 

samples) . 





Taxon 


Net 


Sled 


Entr 


Pontoporeia hoyi 


965.1 


245.0 


1809.4 


Gammarus spp. 


19.6 


36.2 


272.6 


Cranqonyx pseudogracilis 


1.3 


0.3 


215.2 


Hyalella azteca 


0.3 


0.5 


1.1 


Mysis relicta 


50.0 


45.4 


68.8 


Asellus sp. 


0.0 


0.1 


3.7 


Total 


1036.4 


327.5 


2370.8 


n 


476 


140 


542 
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Diel Distribution-- 

Comparison of diel malacostracan densities from net and 
sled tows confirmed the expected nocturnal nature of these 
malacostracans. Of the two methods, diel activity was most 
strongly expressed among net tows, with malacostracan night 
densities being two orders of magnitude greater than day 
densities (Table 19). Net tow diel density differences were 
greatest for M. relicta and P. hoyi , with night densities 
being 336 times and 141 times greater than day densities, 
respectively. As was expected, because sled tows sampled 
very close to the lake bottom in which many of the 
nocturnally active malacostracans take shelter during 
daylight hours, diel ratios were not as extreme as those 
observed in net tows, but nonetheless were indicative of 
diel activity, having an overall night-to-day ratio of 5.2. 

Depth Distribution — 

Very low malacostracan abundances were observed in the 
1- to 6-m depth range when compared with the 9- to 15-m 
depth range sampled by both net and sled tows (Table 20). 

This trend was particularly evident for P. hoyi and 

— 3 ""3 
M. relicta . P. hoyi density averaged 17 x 10 m (sled) 

— 3 ~3 
and 11 x 10 m (net) at depths less than 6 m, but 549 x 

10""^ m""^ (sled) and 153 x lO"*^ m"^ (net) at depths greater 



than 6 m. Similarly, M. relicta densities averaged 4.7 x 
10*^^ m""-^ (sled) and 3.6 x 10~^ m"^ (net) at depths less tl 
6 m, but 100 X 10~^ m"^ (sled) and 78 x 10"^ m""^ (net) at 
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depths greater than 6 m. However, Gammarus spp. and 
C. pseudogracilis densities did not increase with depth, but 
rather were generally similar at depths less than 15 m. In 
addition, once in the water column Gammarus spp., 
C. pseudogracilis , and H. azteca occurred in greatest 
densities near the surface as opposed to mid- and near 
bottom depths in the water column. P. hoyi and M. relicta 
were strongly depth-dependent in the water column, occurring 
most abundantly in the deeper half of the water column. 
This difference may reflect the relative swimming abilities 
of the two groups. P. hoyi and M. relicta are strong 
swimmers capable of a great degree of vertical mobility. 
Remaining taxa are comparatively weaker swimmers and, while 
capable of entering the water column, they may not be strong 
enough swimmers to effectively control their vertical 
position, with subsequent concentration near the surface. 

Seasonal Development and Distribution — 

The seasonal development of P. hoyi based on both net 
and sled tows followed that expected from benthic samples 
collected during April, July, and October, but additionally 
conformed well with the more extensive sampling effort 
conducted near the Cook Plant (unpublished data, GLRD). The 
basic April through September development pattern was as 
follows: primarily spent and gravid females and young-of- 
the-year individuals in April, primarily recently released 
juvenile individuals in May and June (<3 mm), and developing 
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juveniles in July, August, and September (3 to 5 mm) (Table 
21). These data indicated not only the seasonal changes in 
size of P. hoyi individuals, but that net tows, sled tows 
(Table 21), and Ponar grab sampling (Appendix 1) methods 
provided fairly similar estimates of the developmental 
pattern of the population during this time period. 

Seasonally, greatest densities of P. hoyi were found 
during summer months in both net and sled tows (Table 22). 

Peak abundances of P. hoyi occurred during July (3,752 x 

-.0—3 ~3~3 

10 m ) in net tows and during August (625 x 10 m ) in 

sled tows. 

Seasonal distribution of mysids based on net tows 

indicated similar densities occurred April through June 

(approximately 20 to 30 x 10 m ), peak abundance in July 

-3 -3 
(171 x 10 m ), and decreasing density to very low levels 

■"3 ""3 
by September (1 x 10 m ). In contrast, among sled tows 

-3 -3 
lowest mysid abundance was observed in July (4 x 10 m ), 

with August (122 x lO""'^ m""'^) and September (58 x lO""^ m"^) 

densities being highest among the monthly estimates (Table 

22). Whereas the monthly density trends for P. hoyi and 

Gammarus spp. occurring in net and sled tows were similar, 

the lack of monthly comparability for M. relicta is 

suspected to be directly related to the migratory nature of 

mysids in relation to weather, upwelling, and coincidental 

sampling with these events, as well as any non-synchrony of 

sampling effort between sled and net tows. 
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A significant difference between mysids and all 
remaining malacostracans observed in the survey area was 
that the bulk of the mysid population occurred in the 
profundal water of Lake Michigan {>70 m) , thereby 
necessitating migration over a considerable distance in 
order to occur in samples at depths less than 15 m. In 
addition, because mysids are very adept swimmers with very 
good visual perception of light; their ability to avoid mid- 
water sampling devices as well as benthic sampling devices 
is adequate enough to not regard observed densities with the 
same degree of reliability as one would with other benthic 
-macroinvertebrates. Finally, due to the migratory nature of 
mysids, monitoring of seasonal development of mysids caught 
at depths less than 15 m may not be representative of the 
total population. Consequently, little can be stated 
regarding seasonal development of the mysid population. 

Gamma r us spp. exhibited greatest reproductive activity 
during late summer. With occurrence of maximum reproductive 
activity and newly recruited individuals, greatest monthly 
gammarid abundances were observed during August and 
September in net and sled tows (Table 22). 

Entrainment of Malacostracans b y 
Pumping Operations at Unit 3 

General Distribution — 

As was the case in net and sled tows, malacostracan 

entrainment density estimates were dominated by P. hoyi , 

-3 -3 
which comprised 76% (1,809 x 10 m '^) of all entrained 
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malacostracans (Table 18) • Gammarids were the second- and 
third-most numerously entrained taxa, with Gamma r us 

spp. ( Gamma r us fasciatus and Gamma r us pseudol imnaeus ) 

-3-3 . . -3 

averaging 273 x 10 m and C. pseudogracilis 215 x 10 

-3 
m . However, the second-most abundantly occurring 

malacostracan in net and sled tows, M. relicta , was only the 

fourth-most abundantly entrained malacostracan (69 x 10 

As entrained malacostracans originated from Lake 
Michigan, directly comparing entrainment samples (collected 
from January through December) with net tows (collected from 
April through September), to determine how representative 
the former is of the latter, requires examining both 
sampling methods on an April to September basis. In 
addition, it is necessary to restrict net tows to the lower 
half of the water column at a depth of 9 to 15 m as the 
intake structure entrained water from a lake depth of 1 1 m 
at a vertical height above bottom of 1 to 3 m. When these 
restrictions were placed upon the comparisons, extremely 
similar results were obtained, Entrainment of P. hoyi and 
M. relicta averaged 2,347 x 10"^ m*""^ and 101 x 10""^ m""^, 
respectively, from April through September. Average net tow 
density of these species was 2,407 x lO"*^ m""^ and 127 x lO""^ 
m" , respectively. From these comparisons, it was evident 
that entrainment samples were very representative of 
densities in the water column of Lake Michigan in the inner 
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region. However, dissimilar densities were noted among 
gammarids collected by the same two methods, 

A very significant difference between entrainment 
sample and net tow gammarid densities was observed. High 

numbers of Gamma r us spp. (387 x 10 m ) and 

-3 -3 
C. pseudogracilis (308 x 10 m ) were entrained when 

compared with quite sparse numbers in net tows (12 x 10 

m and 0.7 x 10 m , respectively). With respect to 

densities of entrained gammarids, the protective riprap 

covering the bottom area near the intake structures provided 

an ideal environment for establishment of Gamma r us spp. and 

C. pseudogracilis populations. These species occur on rock 

surfaces in several locations throughout the Canadian 

shoreline of the Great Lakes (Barton and Hynes 1976). In 

addition, Gamma r us spp. occur on the riprap at the Cook 

Plant, although Crangonyx has not been observed. Artificial 

substrates placed in the riprap at Campbell were colonized 

by Gamma r us spp., but not by C rangonyx . While Crangonyx was 

not observed on the artificial substrates, its frequent 

occurrence and abundance in entrainment samples strongly 

suggest it has established a large population on the riprap. 

With respect to net tow gammarid densities, the net tow 

sample site was located slightly north of the riprap area 

above sandy substrates. Due to the distance from the riprap 

which supported a large gammarid population, net tows 

collected few gammarids. As the gammarids encountered in 

the study area are best adapted for movement among algae and 
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periphyton growing on rocky substrates, few would be 
expected to occur in an area dominated by shifting sands, 
where movement would be largely limited to burrowing or 
pelagic adaptations common to P. hoyi and M. relicta . Once 
in the water column, both P. hoyi and M. relicta , being 
strong swimmers, burrowers, and migrators, were distributed 
in similar mean abundance over a wide area, irrespective of 
benthic substrate. Subsequently, occurrence of few 
gammarids in net tows, but many in entrainment samples, was 
directly related to behavioral adaptions of the gammarids. 
This suggests that the location of net tows greatly affected 
observed gammarid density, but had no effect on the density 
of either P. hoyi or M. relicta . We conclude that the lack 
of comparability between gammarid densities in entrainment 
samples and net tows was directly attributable to the 
location of net tow sampling in relation to the riprap 
surrounding the intake structures. 

Diel Distribution — 

Entrainment of malacostracans peaked at night and was 
lowest during daylight (Table 19). Dawn and dusk periods 
had similar densities which were intermediate to those of 
the night and day periods. Generally, malacostracan 
entrainment was 16 times greater at night than during the 
day, but ranged from as high as 87 times greater for 
M. relicta and 53 times greater for P. hoyi to as low as two 
times greater for H. azteca . Excepting the former two 
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species, the night-to~day ratio was less than five for all 
remaining taxa. While the absolute ratios differed among 
net tows, sled tows, and entrainment samples, the general 
diel trend was the same. 

Seasonal Distribution and Development — 

Considerable numbers of P. hoyi were entrained during 

June (2,293 x 10^^ m"^), July (8,893 x 10""^ m"-^), and 

~ 3 ""3 
December (5,154 x 10 m ) (Table 23). Maximum peak 

entrainment densities corresponded well with times when one 

might expect greater numbers of P. hoyi to be entrained due 

to recruitment from winter reproductive activity. P. hoyi 

entrained in June and July were the present year class of 

<5-mm individuals recruited into the population during May. 

Those entrained during December were primarily the highly 

mobile, pelagic males present only during winter months, the 

time of maximum reproductive activity (Table 24). Large 

numbers of males were also observed at the Cook Plant in 

December and January (unpublished data, GLRD). Males 

continued to be the dominate reproductive form in 

entrainment samples from December through March at the 

Campbell Plant. Gravid females, which were most evident 

during February, had released their young by April because 

spent females dominated during April. Newly released young 

were first evident in entrainment samples during May. This 

seasonal developnient cycle closely followed that observed 

among entrainment samples collected at the Cook Plant 
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(unpublished data, GLRD) and that observed in net and sled 
tows at the Campbell Plant. 

Mysids were entrained in greatest numbers from May 
through September; density peaked during July (228 x 10'^ 
m ) (Table 23). The majority of entrained mysids were less 
than 6 mm in length, evidently individuals of the current 
year class. Although juvenile mysids were dominant during 
most months, males and females were dominant forms in 
January and February, suggesting maximum reproductive 
activity occurred during winter. 

While the seasonal development of mysids, based on 
individuals collected in entrainment samples at the Campbell 
Plant, was similar to that observed in entrainment samples 
at the Cook Plant (unpublished data, GLRD), seasonal 
distribution of mysids was markedly different. Greatest 
average mysid density at the Cook Plant (1975 to 1978) 
occurred in January during suspected maximum onshore 
reproductive activity, with the majority entrained being 
males. However, mysid entrainment densities were low during 
summer. This general entrainment pattern at the Cook Plant 
was reversed from that at the Campbell Plant. The most 
likely explanation of this difference is that sampling 
effort (during January in particular) at the Campbell Plant 
did not coincide with onshore migration of the mysid 
population to the extent that this event was noted in Cook 
entrainment samples. Entrainment sampling during 1981 at 
the Campbell Plant occurred on 15 and 29 January and on 21 

144 



to 24 December. Quite coincidentally , Cook Plant 
entrainment samples during January 1981 were collected on 
the 12th and 27th. Comparison of the Cook Plant average 
mysid entrainment density estimated from nearly the same 
days with the Campbell Plant average mysid entrainment 

density indicated mysids were entrained in similar densities 

""3 ~3 ~3 

at each plant during January (27 x 10 m and 24 x 10 

"3 
m , respectively) and that onshore-migrating mysids were 

not heavily entrained by either plant during January 1981. 

It is apparent that onshore migration of mature mysids 

occurs from mid-December through January or February. 

Observing an order of magnitude decline in mysid density 

from October to December at >30-m stations in Lake Michigan^ 

Grossnickle and Morgan (1979) speculated that the decline 

may be due to onshore or offshore migration. Onshore mysid 

migration was supported by Cook Plant entrainment data which 

showed that the moderate increase in mysid density in 

December was followed by a large increase in average mysid 

density in January. During February, mysid abundance 

returned to normal levels at the >30-m stations (Grossnickle 

and Morgan 1979). Similarly, entrainment of mysids 

decreased during the February and March time period at the 

Cook Plant to moderate densities which were followed by low 

abundances from April through August. While the correlation 

is in need of more rigorous documentation, there exists the 

potential for demonstrating an early to mid-winter, onshore 

migration and a late winter, offshore migration. 
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Gamma r us spp. were highly abundant during summer and 
fall months (Table 23). Based on presence of males and 
reproductive females, reproduction occurred throughout the 
summer months, with peak recruitment of young during late 

summer. Following the fall months, abundance of Gamma r us 

-3 -3 
spp. decreased notably from approximately 200 x 10 m to 

-3 -3 
less than 10x10 m from December through March. With 

the continual decrease of water temperature and day length 

during fall months, food may be sufficiently depleted to 

induce starvation which, along with predation, may reduce 

the overwintering population to the low levels observed. 



Impact of Entrainment on Malacostracans 
Pontoporeia hoyi — 

Based on the maximum pumping rate at Unit 3 and an 

-3 
average P. hoyi density in entrainment samples of 1.809 m , 

9 - 1 
maximum annual entrainment of P. hoyi was 1.17 x 10 yr 

However, as Unit 3 did not pump at maximum capacity at all 

times during the year, actual entrained P. hoyi on an annual 

basis using the average measured pumping rate (1.115 x 10 

3—1 8—1 2 

m day ) was 7.36 x 10 yr . Assuming a 2-km area in the 

immediate vicinity of the intake structures of Unit 3 as the 

area in Lake Michigan most likely to be impacted by plant 

operations, the total number of P. hoyi available to 

entrainment processes was 1.54 x 10 , based on the inner 

region benthic density at 3 to 15 m during 1981 (3,849 m ). 

Using these density estimates to approximate the potential 
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impact of entrainment , approximately 4.9% of P. hoyi 

2 

occurring m a 2-km area was entrained during 1981. The 

proportion of annually entrained P. hoyi could be as high as 
7.6% were maximal pumping assumed throughout the year. 

Similar calculations utilizing P. hoyi ash-free dry 
weight biomass estimates indicated 518 to 820 kg yr were 

entrained during 1981 at the Campbell Plant. When compared 

2 
with the total benthic biomass available m the 2-km area^ 

annually entrained biomass accounted for 8.4 to 13.3% of 

benthic P. hoyi biomass depending upon the pumping rate 

assumed. 

When compared with similar maximal pumping rates 

8—1 - 1 

estimated at the Cook Plant [1.46 x 10 yr and 297 kg yr 

(ash-free dry weight); unpublished data, GLRD] , the Campbell 

Plant entrained 8.0 times as many P. hoyi as did the Cook 

Plant and 2.8 times the biomass entrained at Cook. The low 

Campbell-to-Cook entrained biomass ratio was at least 

partially due to the method of calculating biomass (Campbell 

calculations used 0.5-mm size classes; whereas, at Cook much 

broader size classes were utilized: <3 mm, 3-5 mm, 5-7 mm, 

and ^7 mm. 

Possibly the inter-plant comparison that best accounts 

for annual pumping rates, entrainment rates, and benthic 

densities of P. hoyi , is the ratio of the benthic area 

required annually to supply the number of P. hoyi annually 

entrained. At the Campbell Plant this area is maximally 

equivalent to 3.04 x 10^ m^ yr"'^ (= 1.17 x 10^ P. hoyi 
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yr /3,849 P. hoyi m ) . At the Cook Plant the similar area 

is maximally equivalent to 1.93 x 10^ m^ yr""^ (= 1.46 x 10^ 

- 1 -2 

P. hoyi yr /757 P. hoyi m ). As the ratio of these two 

numbers compares equally the relative impact of annual 

entrainment on the average benthic density of P. hoyi , 

regardless of entrainment rate, pumping rate, and benthic 

density differences between plants, the relative effect of 

entrainment of P. hoyi was 1.58 times greater at the 

Campbell Plant when compared with the Cook Plant. As 

previously noted, this corresponds maximally to 7.6% of the 

average available P. hoyi population in a 2-km area near 

the Campbell Plant. The equivalent estimate at the Cook 

Plant was 4.8%. 

Given that the effect of entrainment of P. hoyi was 

approximately 50 to 60% greater at the Campbell Plant than 

at the Cook Plant, does this constitute a significant 

difference among relative plant impacts by upsetting either 

the population biology of P. hoyi or the general ecology in 

2 

the 2-km area maximally impacted in Lake Michigan by 

entrainment activities? Neither at the Campbell Plant nor 
at the Cook Plant (unpublished data, GLRD) were any changes 
observed that would indicate an alteration of the population 
had occurred due to entrainment. If it is assumed that all 
individuals passing through the plant were eliminated from 
the population, lake densities were altered. However, as 
benthic abundance of P. hoyi was so high in the area 
considered, we concluded that the loss to the population was 
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negligible. This is especially true if the nearshore 
population is viewed not as a static 2-km area, but rather 
as a churning, swirling, constantly changing dynamism of 
benthic and pelagic environments experiencing both 
immigration and emigration of individuals. During tranquil 
periods the portion of the population near the intake 
structures may experience increased entrainment. 
Conversely, during turbulent periods the effect on any given 
portion of the population may be non-existent due to massive 
movements of the population. On an average day during the 

April to September time period, we calculated there were 

7 
4.61 X 10 P. hoyi m the volume of water overlying the 2- 

2 
km area of concern. The number in the water column 

corresponded to 0.3% of the total available in the water 

mass overlying the 2-km area. Marzolf (1965) estimated the 

percent of the P. hoyi population occurring in the water 

column on 30 different dates in Grand Traverse Bay, Lake 

Michigan at 42 m. While estimates ranged from to 7.4%, 

the average percentage of P. hoyi occurring in the water 

column was 0.96%, suggesting that our estimate was not 

unreasonable. 

As the maximal entrainment rate during this time period 

was 4.15 X 10 P. hoyi on an average day, we determined 9.0% 

(maximal) or 6.0% (actual) of the pelagic population was 

entrained. Overall, we concluded from the 1981 density 

estimates that entrainment maximally affected about 7.6% of 

the benthic population or 9.0% of the pelagic population. 
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but given the ecological dynamics of P. hoyi and the 
physical dynamics of the lake, this effect is diffused over 
a broad spectrum of the population, not continually 
concentrated on a single sub-unit. This perspective, then, 
affects our interpretation of removal of P. hoyi on the 
general ecology in the area. We have assumed the effect of 
entrainment was on a dynamically changing P. hoyi population 
which dominated the benthic habitat at 9 to 15m near the 
Campbell Plant. Under the assumption of a dynamically 
changing population model whereby lake-wide processes unite 
many sub-units of the population along the shoreline of 
eastern Lake Michigan, it is likely that the population is 
continually re-established at rates affected by weather, 
temperature, light, currents, and seasons. Overall, we have 
not observed among these data any trends that would lead us 
to suspect entrainment of P. hoyi has altered the general 
benthic ecology near the plant or the population biology of 
P. hoyi . A possible exception may be the effect entrainment 
exerts on reproductive success during winter months. 

Quite possibly, the time when one might seriously 
consider the effect of entrainment on the population biology 
of P. hoyi is when massive numbers of pelagic males occur 
during peak nearshore reproductive activity in December and 
January. Based on annual entrainment estimates at Campbell 
and Cook, December and January samples were dominated by 
over 90% males. Males entrained during these 2 mo comprised 
17% and 35%, respectively, of the total annually entrained 
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p. hoyi at each plant. If the average number of females 
fertilized by a single male and the average number of eggs 
fertilized by a single encounter were known, the maximum 
loss to the population by entraining a given number of males 
could be calculated, assuming entrained males had not 
fertilized any females and 100% survival of gravid females 
and brooding young until release from the marsupium as 
juveniles. However, without this information, no estimate 
can be generated and the significance of the loss of males 
during the peak reproductive period remains unknown. 

Mysis relicta -- 

Maximum annual density of entrained M. relicta at the 

7 - 1 
Campbell Plant was equivalent to 4.44 x 10 yr or 19.1 kg 

yr (dry weight). Using actual average daily pumping rates 

at Unit 3 from January through December, actual annual 

7 - 1 
density of entrained mysids during 1981 was 2.80 x 10 yr 

or 12.0 kg yr . During the April through September time 

period, we calculated there were 2.39 x 10 mysids in the 

2 

volume of water above the 2-km area during any average day. 

5 
As 1.79 X 10 mysids were maximally entrained daily during 

the same time period, we determined that 7,5% (maximal) or 

5.0% (actual) of the mysids in the water column were 

entrained during the average day. 

Although no benthic density estimates for M. relicta 

were available near the Campbell Plant, Morgan and Beeton 

(1978) estimated an average density of 188 m from their 
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study site in the profundal zone of Lake Michigan near 
Milwaukee, Wisconsin. If this mysid benthic density was 
used to calculate the approximate number of mysids available 

for entrainment at any average time, assuming all were to 

2 

migrate and were entramable, then m a theoretical 2-km 

Q 

area 7.52 x 10 mysids are available. Of this number an 
average 5.9% would be entrained at the Campbell Plant. 
Nearly twice as many mysids were entrained at the Cook Plant 
(8.05 X 10 yr" ) than at the Campbell Plant. The mysids 

entrained at the Cook Plant comprised 10.7% of the 

2 

theoretical 2-km area density at that plant. The greater 

entrainment rate at the Cook Plant was due to the effect of 
December and January sampling which was more strongly 
affected by the mysid density than was the estimate from the 

Campbell Plant. 

. . . 2 

The validity of using a 2-km area as a measure of 

entrainable M. relicta , given the migration distances and 
dynamics of the lake, is certainly questionable. Use of the 
theoretical area could just as easily be increased or 
decreased, with subsequent effects on percentages entrained, 
not just for mysids, but for amphipods as well. Its use 
serves only one purpose: to add perspective to the magnitude 
of the numbers or organisms entrained and those potentially 
entrainable. No absolute is implied by the percentages 
generated. However, it is interesting to note the extremely 
large numbers of mysids that are available for entrainment 
assuming 100% onshore migration. Without knowledge of 
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approximate benthic mysid densities at a site, the relative 
effect of entrainment at sites can not be fully evaluated. 
Current benthic density estimates for the mysid population 
are poorly known. In addition, very little is known about 
the degree of onshore migration by the mysid population. 
Consequently, the effect of entrainment on the population of 
M. relicta is even more poorly understood. 

Gamma r us spp. and Cranqonyx pseudoqracilis — 

While annual entrainment could be determined, 
estimating the effect of entrainment on the populations of 
G. fasciatus , G. pseudol imnaeus , and C. pseudoqracilis was 
specific to the month of October as the only estimate of 
benthic density was from artificial substrates retrieved 
during October. Annual entrainment of Gamma r us spp. was 

estimated to have a maximum of 1.76 x 10 yr , while that 

8 — 1 
of C. pseudoqracilis was 1.39 xlO yr . However, based on 

the actual average daily pumping rate at Unit 3, actual 

8 "- 1 
annual entrainment was 1.11 x 10 yr for Gamma r us spp. and 

7 - 1 
8.76 X 10 yr for C. pseudoqracilis . No biomass estimate 

was determined for C. pseudoqracilis , but maximal and actual 

annual ash-free dry weight biomass estimates for Gamma r us 

spp. were 50.8 kg yr" and 32.1 kg yr" , respectively. As a 

good annual benthic gammarid density estimate was not 

available from the riprap from which the gammarids 

originated, the effect of entrainment was based on the 
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October entrainment rate (0.2005 m ) and the average 

-2 

October benthic density of 139 Gamma r us spp. m 

2 
In the riprap area, which is equivalent to 52,400 m , 

7.28 X 10 Gammarus spp. were potentially entrainable during 

October. Using the maximal pumping rate at Unit 3 and the 

5 
October entrainment density, an average 3.54 x 10 Gammarus 

spp. were entrained daily at the Campbell Plant in October. 

This estimate corresponded to 4.9% of the available riprap 

population. Using the actual daily pumping rate at Unit 3 

in October (94.8 x 10^ m^ day"^), 1.90 x 10^ Gammarus 

spp. were entrained, or 2.6% of the total available riprap 

population. 

-3 
Although C. pseudogracilis averaged 0.1306 m m 

October entrainment samples, none were observed on 

artificial substrates. The lack of Crangonyx on the 

artificial substrates and its abundance in entrainment 

samples suggested the artificial substrates may have been a 

poor substrate for Crangonyx to colonize, possibly due to 

physical location on the riprap or inability to support 

suitable food sources. Nevertheless, the large numbers of 

C. pseudogracilis entrained on an annual basis (8.76 x 10 

yr ) and its very poor representation in net and sled tow 

samples make it highly unlikely that entrained Crangonyx 

originate from some other portion of the survey area outside 

the riprap area. Rather, a population nearly equivalent to 

that of Gammarus spp. must reside on the riprap which was 

most certainly the source of entrained C. pseudogracilis . 
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We concluded that artificial substrates did not 
adequately estimate the population densities of the 
gammarids encountered on the riprap. This conclusion was 
reasonable for Gamma r us spp. because within a short time, 
the daily entrainment rate for Gamma r us spp. would deplete 
the estimated riprap population density of Gamma r us spp. to 
near zero at a declining exponential rate. Despite high 
entrainment densities, absence of C. pseudogracilis on 
artificial substrates suggested the artificial substrates 
did not adequately estimate the Cranqonyx population 
occurring on the riprap. 

Hyalella azteca and Asellus sp. — 

Occurrence of H. azteca was limited to only a few 
individuals in net tows, sled tows, and entrainment samples. 
No specimens of H. azteca were found in riprap or Ponar grab 
samples. Averaging only 0.0011 m" in entrainment samples, 
the maximal annual entrainment estimate was 7.09 x 10^ yr~\ 
The actual entrainment estimate of Hyalella was 4.48 x 10^ 
yr . From these data it is apparent that conditions in the 
survey area did not favor establishment of a dominant 
population on the riprap, at least in comparison with the 
gammarid population. A similar situation was observed at 
the Cook Plant where H. azteca averaged only 0.0008 m""'^ in 
entrainment samples (unpublished data, GLRD). 

Abundance of entrained Asellus sp. was 0.0037 m""^, 
which was equivalent to a maximal annual entrained density 
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of 2.38 X 10 individuals. The actual annual entrained 
density for asellids was 1.51 x 10 individuals. 
Entrainment of Asellus at the Cook Plant averaged slightly 
lower (0.0017 m ) than the present estimate at the Campbell 
Plant. Asellids did not occur in net tows and were 
extremely rare in sled tows (0.0001 m ) at the Campbell 
Plant. 

As asellids are cryptic and remain at or very near the 
substrate on which they occur, poor representation in net, 
sled, and entrainment samples was not unexpected. The 
presence of asellids in Ponar grabs, which sampled sandy 
substrates, was limited to only rare occurrences in the 
inner region near the riprap. 

Asellus sp. was the most abundant malacostracan 
colonizing the artificial substrates (313 m ). Average 

asellid abundance on the riprap during October was 

7 
equivalent to 1.64 x 10 individuals. When compared with 

6 3 
the number maximally entrained in October [(1.766 x 10 m 

day"'b(3.079 x 10~^ individuals m"-^) = 5.44 x 10^ 

individuals daily], maximal daily entrainment represented 

only 0.003% of the asellid population occurring on the 

riprap. The very low numbers entrained underscore the 

cryptic nature of asellids and that they are not as pelagic 

as other malacostracans encountered at the Campbell Plant, 

thereby reducing their susceptibility to intense 

entrainment. 
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SUMMARY 



Macrobenthic invertebrates were collected from Ponar 
grab samples taken in the 3- to 15-m depth regime during 
April, July, and October 1978 through 1981 near the 
J,H. Campbell Plant, located along the eastern shoreline of 
Lake Michigan. In addition, from April through September 
1981, malacostracans were enumerated in net and sled tow 
samples collected from vertical depth strata (net) and near 
bottom (sled) at stations ranging from 1 to 15 m in depth. 
Malacostracans occurring in entrainment samples were also 
enumerated from January through December 1981. Finally, 
density was determined for macroinvertebrates collected from 
artificial substrates lodged in the riprap during July 1981 
and retrieved during October 1981. These data provided the 
information from which we were able to determine the effect 
of discharged heated effluent and entrainment on various 
major constituents of the macroinvertebrate community common 
to the Lake Michigan benthic environment near the Campbell 
Plant. 

An analysis of variance (ANOVA) was performed on each 
of the nine major macroinvertebrate constituents of the 
benthic community, plus total benthos. In order of 
decreasing numerical importance in the survey area, the nine 
most numerous macroinvertebrates were: Pontoporeia hoyi , 
Chironomidae, Naididae, Tubificidae, Pisidium , Turbellaria, 
Stylodrilus herinqianus , Enchytraeidae, and Gastropoda. 
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Effect due to discharge of heated effluent was determined by 
comparing R (least detectable true difference; calculated 
using ANOVA mean square error) with R' (actual abundance 
change ratio; calculated from averaged regional log 
densities before and after plant operation) (see METHODS for 
detail). Based on this comparison for each taxon, we 
determined no detectable effect due to heated effluent 
discharge was attributable to density fluctuations observed 
among the 10 groups tested. Therefore, we have concluded 
that operation of the Campbell Plant through 1981, with 
subsequent offshore discharge of heated effluent, has had 
either no effect or no detectable effect on those individual 
portions of the macrobenthic community comprising over 99% 
of the total benthic density or on the total benthic 
density, regardless of regional or annual density 
differences noted from 1978 through 1981. 

Numerically, the malacostracan most often entrained was 
P. hoyi ; followed in decreasing numerical order by Gamma r us 
spp. (G. fasciatus and G. pseudolimnaeus ) , Cranqonyx 
pseudogracilis , Mysis relicta , Asellus sp., and Hyalella 
azteca . On an annual basis, maximal and actual (as 

determined from maximum and actual daily pumping rates at 

9 
the Campbell Plant) entrainment of P. hoyi was 1.17 x 10 

yr"^ (820 kg yr"\ ash-free dry weight) and 7.36 x 10^ yr"^ 

(518 kg yr~ , ash-free dry weight), respectively. Similar 

7-1 1 

values for M. relicta were 4.44 x 10 yr (19.1 kg yr , dry 

weight) (maximal) or 2.80 x 10^ yr""^ (12.0 kg yr"\ dry 
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weight) (actual). When comparison between entrainment 
sample and net tow sample densities was made from similar 
depths and months during which samples were collected using 
both methods, very similar densities between the two methods 
for both taxa were observed, indicating entrainment samples 
were very representative of the pelagic, nearshore 
densities. This trend provided assurance that annual 
densities of entrained organisms were relatively unbiased 
and reasonably accurate estimates of the densities of lake 
organisms potentially entrainable. 

Likely the most important comparison using entrainment 

data was that between number annually entrained and number 

2 

potentially available in a 2-km area. For P. hoyi the 

number potentially entrainable was equivalent to 1.54 xlO 
individuals. Based on the ratio noted above, at most 7.6% 
of the P. hoyi population was entrained during any average 
time period. As the population is very mobile, the effect 
of entrainment on P. hoyi was considered negligible, with 
the possible exception of an unknown effect due to heavy 
entrainment of males during winter reproductive activity in 
December and January. 

The effect of entrainment on M. relicta was also 
negligible due to the high number of mysids potentially 
available from the vast profundal area of Lake Michigan 
immediately offshore from the Campbell Plant. Estimating an 
average mysid density of only 188 m""^ (Morgan and Beeton 
1978) from the offshore in a 2-kin area and assuming all are 
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potentially entrainable via onshore migration, maximally 
only 5.9% of these mysids would be entrained annually. 
However, as with P. hoyi , little is known about the effect 
of entraining large numbers of males during winter 
reproductive activity. 

Remaining malacostracan taxa in samples collected near 
the Campbell Plant originated nearly exclusively on the 
riprap surrounding intake and discharge structures. 
Although entrained Gamma r us spp. and C. pseudogracilis 
occurred in very similar densities, the latter species was 
not found on artificial substrates. High densities of 
C. pseudogracilis in entrainment samples and low to no 
abundance in all samples from remaining collection methods 
left little doubt that their origin was the riprap area. Of 
the two remaining malacostracans encountered, neither 
Asellus sp. nor H. azteca were entrained in high abundance. 
While the latter species did not occur on artificial 
substrates, the former was the most numerous malacostracan 
on these same substrates. The effect of entrainment on 
asellids was very minimal (0.003% of the riprap population) 
and was likely due to the clambering nature of Asellus 
sp. as opposed to the swimming nature of other 
malacostracans encountered; susceptibility of asellid to 
entrainment was thus reduced substantially. 
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